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Abstract. Consider the quintic complex Ginzburg-Landau equation 
ut ={m + ia^j,o)uxx — (in + iafj.i)u 

+ (1 + iafi2)\u\'^u — (1 + ia^3)\u\'^u, a; € R. 

The parameter m is close to 3/16 so that for a = 0, this equation possesses 
a unstable pulse-like solution. For |a| small the equation possesses pulse- 
like solutions of the form e*"'e"^'^V(a;), with r a positive function decreasing 
exponentially at infinity and (p asymptotic to — C|a;| -|- D at infinity. 

These solutions are linearly unstable for \a\ < olc\ if a certain rational 
function of and /^s is strictly positive and not too large, they become stable 
for |a| > «c: when the initial data is a pulse plus a small perturbation, its 
limit for large times is the same pulse, possibly translated in space and in 
phase. 

This article gives a rigorous proof of a conjecture of Thual and Fauve ; it 
relies on a very detailed asymptotic analysis of the eigenvalues of the linearized 
operator, depending on the parameters 3/16 — m and ^ij. 



1. Introduction 

In Thual and Fauve proposed a model of localized structures generated 
by subcritical instabilities; in their article, they mentioned several examples of 
such localized structures in systems far from equilibrium: local regions of turbulent 
motion surrounded by laminar flow as in |^ chapter 19, spatially localized traveling 
waves at convection onset in binary fluid mixtures as in |Q or a Faraday 
experiment in a narrow annular dish as in [ p9[ . 

The equation proposed by Thual and Fauve is the quintic Ginzburg-Landau 
equation 

(1.1) — = mo^^ + miM + m2|upu -I- TO3|u|''u, 

where the rrii are complex coefficients with Jftmo > 0. 

Thual and Fauve assumed 5Rm3 < 0, in order to stabilize large amplitudes, 
and SRtoi < 0, so that the zero solution should be stable. Choosing 5ftTO2 > 
and adequate relations on the other parameters of the equation ensured that there 
would exist non zero stable homogeneous solutions of the form 7.e*'^*+*''^. 



They integrated numerically (1.1) and found thus a solution of (1.1) of the form 



(1.2) u(a;,t) = e*"*r(a;)e*'^(^\ 

where r takes positive values and decays exponentially at infinity, while the phase (j) 
is asymptotic to —C|a;| -I- D at = oo. They obtained numerically standing wave 



solutions of (L3), which they found to be experimentally stable for large enough 



absolute values of the imaginary parts of toi, 7712 and TO3. 
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Let US first simplify (LI) tlianks to some scale consideration: by clianging the 
units of t, u and x, we can see that it is possible to choose 5Rto2 = 1, 5Rm3 = — 1 
and m = Sftwo = — 3?toi. Thanks to the equivariance, we can also take toq to 
be real, with very little loss of generality. For reasons which will be justified in the 
course of the article, our results hold when m is slightly less than 3/16; let it be 
said only at this point that for m — 3/16, and all the m^'s real, there exist four real 
heteroclinic solutions of 

—mr" + mr — + — 0, 

which are distinct up to translation. We let f be one of these heteroclinic solutions 
which takes the value y/3/2 at ~oo and at +oo; the three other ones are obtained 
by mirroring f, a; or both. 

For m slightly less than 3/16, there are two homoclinic solutions with a very large 
"shelf", i.e. a region where the solution is very close to ±3/4; it is the presence of 
this large shelf which makes this choice of parameters interesting. 

For real values of m^, pulses must be linearly unstable. Let us sketch the ar- 
gument which proves this statement: if u is a pulse- like solution according to the 



above definition, it will be proved at lemma 2.1 that u can be taken real, even and 



positive; then u solves the following differential equation: 
(1.3) -muxx + mu - + ^ 0, xeR. 

Moreover, the linearized operator at u is the unbounded operator in i^(M) given 
by 



(1.4) D{A)=H'\R), As^ -ms" + ms-3u^s + 5u'^^ 



s. 



The operator A can also be seen as a Schrodinger operator in M, with potential 
m— 3m^-|-5u*; in particular, its essential spectrum lies above m. Differentiating ( |l.3| ) 
with respect to x, we find 

Au' = 0, 



which expresses the translation invariance of (1.3). Thus u' is an eigenf unction 
relative to the eigenvalue 0; since u' changes sign, the maximum principle implies 
that cannot be the lower bound of the spectrum of A. We denote by A < this 
lower bound. 

The eigenmode pertaining to A is the shrinking-swelling mode: under most per- 
turbations, the pulse shrinks to or swells to infinity. 

In order to prove our main results, we introduce a few notations: 

16m 

(1.5) ^=i__, 

1 4 

(1.6) ^ = ll^^- 

4 v 

In the first part of this introduction, we restrict ourselves to a simplified case 
of (|li|): 

(1.7) Ut — ^mu + mUxx + (1 + ia)|upM — 

where m and a are real parameters; 

In this article, we prove the following existence theorem: 



Theorem (Theorem 4.5). For all p > 0, for all v small enough, there exists C > 



such that for all a satisfying 



\a\<a„, = ]^^{l-CL-P) 
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there exists a pulse solution of (L7), i.e. a solution of the form 

uix,t) = e*'^*r(a;)e''^(^\ 

where r is a positive function which decays exponentially at infinity, and 4>{x) is 
asymptotic at infinity to —C\x\ + D, with C a positive constant and D a real con- 
stant. 

Then, we prove the stabiUzation property conjectured by Thual and Fauve: 

Proposition (Proposition 5.8). Let u he the solution defined at theorem 4-5. Let 
\o^\ < ctm be the above defined number; there exists a number ac whose asymptotic 
is 



1 



48L^ 



such that the solution u is stable iff olc < \oi\ < am, and unstable if \a\ < ac- 

More precisely, if \c(\ < ac, the spectrum of the linearized operator at u con- 
tains exactly one negative eigenvalue; when a = ac, the eigenvalue is of algebraic 
multiplicity 3 and geometric multiplicity 2, with a non trivial Jordan block of di- 
mension 2; when ac < \a\ < am, the linearized operator at u has the semisimple 
double eigenvalue 0, and the remainder of the spectrum is bounded away from the 
imaginary axis. 

The general case i.e. 

(1.8) Ut = rnuxx — (m + iafii)u + (1 + ia^2) u — (1 + ia^a) \u\^ u. 

is treated in section with very few analytical details; under the condition 



X(a) = 



Ai2 



9/^3 



3^ 
16 



9m3 
16 



1 



{2^12 ~ 15Ai3/8)' 



>0, 



skew stabilization also takes place, as is shown in Proposition |6.l| . 

Let us give a very rough idea of the reason for existence, and of the skew sta- 
bilization mechanism. In both cases, we will restrict ourselves to the simplified 
equation (1.7). 

For the existence, we start from so-called kinks: they are solutions of (1.7) of 
the form 

(1.9) u{x,t) = e"^*K{x-ct), 
with 

(1.10) K{x) ^r{x)e'^'=. 

to make things precise, we demand that r take p ositi ve val ues, i ncreasi ng from 
r{-oo) = to f = r(+oo) > 0. If we substitute (1^) and (lAX.) into (L7), we 
obtain the equation 



m(r" + 2ikr' — k'^r) -f mr — (1 4- ia)r'^ + r^ 



0. 



(1.11) ibjr ~ c{r' + ikr) 
Define uj — uj ~ kc; the vanishing of the imaginary part of (1.11) implies that 



(1.12) 



- ar^)r 
~2k ' 



Since we assumed that r is increasing from to r and is strictly positive, rela- 
tion (1.12) implies the sign conditions: 

(1.13) Cj/k>Q, a/k>0. 

If we differentiate this relation with respect to x, we obtain 

„ {uj — ar'^){uj ~ 3ar'^)r 



(1.14) 



4fc2 
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Substituting the expressions ( 1.12 ) and ( 1.14 ) in the real part of ( 1.11 ), we obtain 



a polynomial of degree 5 in r; since r is assumed to be different from 0, we obtain 
through straightforward algebra 

^ - — , 2Cj + ck=—, Lu'^ -aCj + a^(m + ia^/A) = Q. 
4to 2 

The second degree equation in uj has real roots if and only if Aa^ — v < 1/3, which 
we assume from now on. These roots are given by 

(1.15) Cj^j{2±^/1 + 2,{v- 4a2 

We see also that 



(1.16) f=V^^- 

If we linearize the equation for v = u{x — ct) exp(— z/c(a; — ct) — iujt^ around f, the 
condition for stability of 0-wave number modes is 



2f2 _ 4^4 ^ 



in (1.15) 



Co mpar ing this relation with ( 1.16 ), we can see that we have to choose the + sign 



Hence, the velocity c is given by 

\/3(4a2 _ 



1 + Vl + 3(i^-4a2) 



where we have used the sign condition (1.13). 

If 4a^ > V, the velocity c is positive, and the zero state gains over the non zero 
state; we will say that this is the inflow situation; on the contrary, if Aa^ < v, 
the non zero state gains over the zero state; we will say that we have an outflow 
situation. 

Assume from now on that a and v are small. If there existed a pulse- like solution, 
it could probably be approximated by a combination of a kink K{x + L — ct)e^'^* 
centered at —L and an antikink K{L~x^ct) centered at L, provided that we know 
how to glue their phases together. In the outflow case, two competing effects take 
place: on one hand, the choice of parameters tends to produce an expanding pulse; 
on the other hand, an attraction effect between the walls is expected as in [||, |l^ 
and jl^ . It is natural to expect that this attraction effect should be exponentially 
small with L. Thus, it is reasonable to conjecture that the evolution of the pulse 
will be given by 

L = -c-Cie-^^^. 
The pulse will be in an equilibrium if the two competing effects balance, i.e. 

This analysis assumes an almost scalar pulse, so that it is easy to glue together 
the phase of the kink and of the antikink. Of course, the pulse obtained by this 
argument is not stable: any amount of swelling or shrinking of L destabilizes it. 

Let us consider the inflow case: Aa^ > v. Now, nothing scalar can stop the 
pulse from collapsing. However, Malomed and Nepomnyashchy argue in ]3C| ] that 
the pulse does not collapse because of phase incompatibility: a minimum distance 
is necessary to match the phases of the kink and of the antikink. By formal asymp- 
totics, they obtain a half length of the pulse given by 

in a very small range ^ Aa^ — <C a^. 
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Figure 1. The length L of the pulse as a function of a.. 



However, the second author of the present article tends to believe that ( l.lq ) 
holds for a much larger range of a. The article |30| contains no analysis of the 
stability of the pulse obtained; however the claim of stability (in a mathematical 
sense) seems quite reasonable. 

With our results and the results of Malomed and Nepomnyashchy, we can plot 
a graph of L as a function of a, and we obtain Fig. |^. 

The region around — v seems difficult and interesting; jumping somewhat 
too fast to conclusions, it would be nice to believe that these two branches join to 
form a single branch; proving this might mean considerable effort. 

Let us now give an idea of the mechanism of skew stabilization. We recall 
that we assumed that the shelf in the solution is very large, i.e. v is very small. 
The definition of L is such that for a = 0, the width of the shelf of the pulse is 
approximately 2L. We define a function 



F(u^ uj, a, v) = iuju - 



mu 



+ mu — (1 + iQ;)|u| u + \u\ 



A function ue*'^* is a solution of (1.7) iff u and uj are such that F{u,LO,a,v) van- 
ishes. Moreover, it is linearly stable if DuF{u,uj,a,v) has its spectrum included 
in the right-hand side complex half-plane; moreover, we require for stability that 
should be a semi-simple eigenvalue of finite multiplicity while the remainder of the 
spectrum is bounded away from the imaginary axis. 

The symmetries of the problem imply that there are two eigenfunctions of D.^F , 
relative to the eigenvalue 0: t ^ iu and u'; there is also the eigenfunction w 
relative to the shrinking-swelling mode, and the aim of the game is to show that 
the corresponding eigenvalue A crosses for some appropriate value of a. Let us 
define 

F{u, T, L, (3) ^ F{u, e~^^(3T, €"'^^(3, e"*^) . 

At this point, we will use a measure of cheating in order to explain what is going 
on: if u and t can be seen as smooth functions of L, we differentiate with respect 
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to L, and wc find 



(1.19) 



D.^F{u{-,L, (3),t{L, /?), L, (3) = e-^"" + 2ze-2^/3(™ - \u\^u) + 

r 2L;3_.. I ,|2„, , I ,|4„, •„-2L, 



Since we expect the interaction between the rotation mode and the shrinking- 
swelling mode to be the reason for the skew stabiUzation, we let {w,t} be a dual 
basis to the basis i.e. 

t tdx= w wdx = l, It w dx = w tdx = 0, 



y and t are even, i"^ DuF{u{L, P),t{L, 0),L, f3) ~ 0, 
w^DuFiuiL,f3),T{L,f3),L,l3) = w^X{L,(3). 
We multiply ( |1.19| ) on the left by , we integrate, and we find that 



X{L,f3) J w^^dx = -2e-~^^l3 J i 



w'^i\u\'^udx 



y"u;^[-|wpu+ |u|4u_ie"2L^|w|2w] dx. 



Since a is small, the problem is very close to being self-adjoint, it is not unreasonable 
to take w as an approximation of w; moreover, in the approximation of the large 
shelf, u is very close to r{\x\ — L); therefore, w and du/dL can be taken very close 
to — f'(|a;| — L)sgna;. These considerations imply that 



w'^ —— dx ^ 2 f If'l^ = — . 



dL 



L 



We see also that 



lupu + lul^u] dx^2 I f'f^ - f'f^ dx . 

64 



L 



Thus, we obtain the following "equivalent" for X{L,P): 

X{L,f3) 2e-^^f3 / s'^i\ufudx. 

If we are able to calculate with enough precision the integral in the right hand 
side of the above equation, and if it turns out to be negative, we will hope for a 
stabilization effect. 

However, we have cheated too much for this argument to go through a rigorous 
analysis; the main source of inexactitude comes from the assumption that for a 
small, we can find an even solution u of F{u, r, L, (3) which is close to the real even 
solution of F{u, 0, L, 0) — 0; indeed, u can be found close to the real even solution 
of F{u, 0, L + 2/, 0) = 0, where y is related to a by the relation 

2/ = -^ln(4/3-2) =-iln(4e2^a-2). 

In the course of the proof of skew stabilization, we will discover the following spatial 
scales in the problem: we have already 1 and L; the third scale is InL, which is the 
order of magnitude of the critical y for which the skew stabilization occurs. 

We can now explain the organization of the paper; since we are in an almost 
scalar situation, we let a — -Je and we define 

G'(^, 77, r, e, v) = 3?F(^ + i\ler\, \feT, ^/e, v) + + i\/eri, y/sT, -/e, v). 
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In our continuation analysis, we need a starting point {u,uj) such that 
(1.20) F{u,u;,0,iy) ^0. 



In section we prove that a solution u of (1.20) exists and decays exponentially 
at infinity only if m belongs to the interval (0, 3/16). This solution is, up to phase 
and space translations, the unique positive even solution of 

— mr" + mr — + = 0, 



as is proved at Lemma 2.1 



Then, we proceed to study precisely this r — r{x,i'), which has an explicit 
expression; however, we use mostly the asymptotic for r^x^v), when j/ is small 
(Le mma 2.2). Thus, we see that the shelf of r has indeed a half-length of L defined 
by (L6). We study the linearized operator A around r; it is an unbounded self- 
adjoint operator in L^(R) given by 

D{A) = H'^{R), Au = -mu" + (to - Sr^ -|- hr^)u. 

We prove that the spectrum of A contains a group of eigenvalues {A, 0}, and that 

is the lowest eigenvalue of A\ moreover, this group of eigenvalues is bounded away 
from the remainder of the spectrum, uniformly in v. The eigenvalue corresponds 
to the translation mode, with eigenfunction r', and the eigenvalue A corresponds to 
the shrinking-swelling mode, with eigenfunction s. This analysis is made possible 
by the following fact: if 

dr{x,v) d 

cr(x,j/) = — = —r(x,i^}, 

ov oL 

then, with an appropriate normalization of s. 



as is proved at Theorem 2.8. 

The existence will be proved using the scaled equation 

G(f,?7,T,e,j/) = 0. 

The scaled equation is more interesting from the point of view of continuation, 
because for a = 0, cu = 0, and it is easy to check that 

^F{r,0,0,iy) = 0, 

which does not provide any information. However, 

5RG(e,r?,T,0,j.) =5RF(C,0,Q,j/) 

so that r satisfies 

3fiG(r,77,r,0,i/) = 0, 

and 77 and r are yet undetermined. The second equation is 

(1.21) 3G((r, r], r, 0, v) = rr — — mrj" + mrj — r^rj + r'^rj = 0, 

and it is studied in details in Section |^. Let 9 and q denote the values of ij and r 
which satisfy ( 1.21 ). It is natural to study the operator B in i^(R) defined by 

D{B)^H^{M.), Bu ^ ~mu" + {m ~ + r'^)u. 

We can see immediately that is the lowest eigenvalue of B and the correspond- 
ing eigenvector is r; this is not surprising, since it is the analytical translation 
of the equivariance of (1.7). In other words, ir is the phase rotation mode. 
Therefore, (1.21) will have a solution if and only if r is orthogonal to — Or; this 
determines 9] if we impose that q be orthogonal to r, it is uniquely determined. 
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For later purposes, we need an asymptotic on the second eigenvalue ^2 oi B; it 
is proved at Theorem 3.4 that ^2 ^ C/L"^, where C is a positive constant. 



Section ^ is devoted to the existence proof. Preliminary computations showed 
that continuation is not good enough to obtain a satisfactory range of existence; 
what is needed is an ansatz for the pulse; it is obtained by taking — ve~^'^ , 
where y is some positive number bounded by L^] the corresponding r, 6 and q 
are denoted by , 9^ and . Now, has to be determined: this is a version of 
the Lyapunov-Schmidt method of bifurcation theory. Our choice is to require that 
^G{r^ ,0^ ,v) is orthogonal to s^, the shrinking-swelling mode corresponding 



to v" . An asymptotic for e is given by 



4' 



e 

where 

K^{V-V^)/{1-V^). 

Let be the vector of components (r^, q^ ,6^ ,e^). The idea is to observe that this 
is almost a solution of G{U^ , = 0. Existence is proved using a version of the 
implicit function theorem with estimates, proved in Section ^ . The pulse obtained 
this way is denoted by u. 

In other words, we approximate the pulse at a = and by the pulse at 
a = and . 

Let us denote by 0(1) any quantity bounded by a finite power of L. 

The main result of this article, i.e. the proof of stabilization (section |^) uses the 
details of the proof of existence. If (1 — k)!) is the differential of F with respect to 
u at (M,a;,a, :^), then. 

Here, 



A' 



A' 




where (resp. B^) \s A (resp. B) at instead of i/, and C are 2x2 matrix of 
multiplication operators C^, 1 < i,j < 2 such that 



\\Bi2\\l^ + IIS21IL00 + IICnILo. + IIC22IL00 = 0{l). 

The idea is to consider the restriction of T) to the generalized eigenspace cor- 
responding to the eigenvalues of B which are close to zero. This eigenspace is of 
dimension 3; a basis of it is {s,iu,u'}, where iu spans the phase rotation mode, 
u' spans the space translation mode, and s corresponds to the shrinking-swelling 
mode. In this basis, the matrix of B is given by 

^Mii 0^ 
M21 
Oy 

The sign of Mn determines the stability of the pulse: if Mn > 0, the pulse is 
stable (up to space and phase translation); if Mn < 0, the pulse is unstable. Thus 
we have to give an asymptotic for Mn. For this purpose, we embed the operator 
C into an holomorphic family "Die) of operators depending on c, and we prove 
estimates using the strong properties of such families. In particular, we give a 
precise description of the expansion of a basis of eigenfunctions relative to the very 
small eigenvalues of 2?(c), and of the dual basis, and we validate these expansions. 
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With the residue theorem, we are able to describe Mu with sufficient precision, 
and the symmetries of the problem lead us to an almost completely explicit value 
for it (Lemma 5.4, Theorem 5.5 and Theorem 5.6). We conclude this asymptotic 
analysis at Proposition 5.S. 

Thanks to a result of Henry ||l^, Chapter 5, Exercise 6, the linearized stability 
implies the following non linear stability result: take an initial condition for (1.7) 
which is equal to a pulse plus a small perturbation ; then, if \a\ > ac, the asymptotic 
state of the solution of ( |l.7[ ) is a pulse possibly translated in space and in phase. 

In section ^ we give the analogous asymptotic for the case when the /ij's do not 
vanish; the appendix (section |^) gives an implicit function theorem with estimates; 
this theorem is the key to the existence result; in other words: our existence result 
is based on an ansatz: if the ansatz is good enough, then it is indeed a good 
approximation of the solution. When small parameters are involved, a correct 
argument deserves a proof. 

A rather curious fact is that the number 7r^/6 appears in the calculation of the 
expansion of Mu; in partial differential equations, it is usually related to a trace, 
but we have been unable to uncover such an origin; therefore, its presence may be 
a coincidence. 

Some of these results were announced in |Q which contains a number of errors. 
A preprint |^ was circulated but never published as an article; the present article 
contains for the first time the approximate explanation of the skew-stabilization of 
the pulse and also the case of general coefficients as in ( |l.8| ). 

There is considerable interest in the Ginzburg-Landau models; scanning the liter- 
ature, one can find for example [^6| which lists fluctuations in lasers, order-disorder 
transitions, population dynamics and ordering in uniaxial ferromagnetic films as 
domains where Ginzburg-Landau of the third degree has been used as a model. 
W. Eckhaus |^ states that Ginzburg-Landau of the third degree is "universal" for 
modulation equations, which is another way of saying that it behaves as a normal 
form. 

The article describes a large number of solutions of (1.1), perturbation 
expansions for large values of 9toj and gives conjectures on the behavior of the 
solutions of (1.1) in different regions of the parameter space. 

Ginzburg-Landau of the fifth order is much less generic than the third order 
Ginzburg-Landau. Its main merit is that it allows for subcritical bifurcation of the 
constant amplitude solutions. 

Thual and Fauve explained in | pl| the phenomenon they obtained in terms of the 
general picture of subcritical bifurcation, and also as a perturbation with respect 
to a nonlinear Schrodinger equation. 

Sh ortl y later, Malomed and Nepomnyashchy considered the same equa- 
tion (IT) and explained the existence of a pulse by formal asymptotics. A careful 
examination of their results shows that they worked in a different range of param- 
eters from ours. 

Hakim, Jakobsen and Pomeau have given a general idea of the bifurcation 
picture in a situation which is close to the present one; however, it is difficult to 
compare the situations, since their statements are not described with complete pre- 
cision. One of their statements is the subcritical character of the bifurcation. While 
the bifurcation of space independent solutions is clearly subcritical, the bifurcation 
of the pulse is subcritical only by the fact that initially the solution is unstable, and 
it is stabilized along the solution curve; however, the typical picture of subcriticality 
as in Fig.l o f has not been found in the present article. 

Kapitula [^2| gave a general theory for the existence of heteroclinic traveling wave 
solutions of the quintic Ginzburg-Landau equation with convective terms. Kapitula 
basically studies the persistence under perturbations of heteroclinic orbits close to 
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orbits which can be obtained ahiiost exphcitly when all coefficients are real. The 
proof relies on a very precise study of the perturbed invariant manifolds for a flow 
associated to the system. Poss ibly, his methods could be adapted to give existence 
of homoclinic solutions of ( |l.3| ). 

According to a very striking phrase by Yves Meyer [ pT| , the present work is like 
the success of the lock breaker: he/she has to use many pick locks, and the place 
looks messy; however, it is expected that once the door is open, and the lock can be 
dismantled and studied in details, someone will be able to devise a nice key which 
will open it in a single move. 

Kapitula used an entirely different set of techniques to devise a nice key for 
related existence questions, but a nice key to stability does not yet exist. 

This article makes available a number of pick locks, to be stored in the tool box 
of the mathematical lock picker. It is more in the spirit of the SLEP method of 
than in the spirit of the large literature on the analysis of the stability of traveling 
waves: solutions of the quintic Ginzburg-Landau equation have been analyzed in 
]23[ |, psj , | p7[ |, p^ , where the main difficulty is the bifurcation from the essential 
spectrum; the foundational work on the stability of traveling and standing wave 
solutions of semilinear parabolic equations is related to the Evans function; see 
in particular [|l9), and which were followed by a considerable literature, 
including in particular [po| , |Q, |2^, The analysis of perturbation of 

periodic states has been taken up by the Evans function method in [Q and by 
modulation equation methods in |40(| . Solutions with several fronts or bumps have 
been constructed in @, [|, ||, |f, g and H. 

It is possible that the methods of this article are close to those of which 
however does not have the symmetry of the problem considered in the present 
article. 

We would like to thank Stephan Fauvc for introducing us to this problem, which 
revealed itself as infinitely more complicated than what we would have expected. 
The second author is glad to thank R.L. Pego, B. Malomed and A. Nepomnyashchy 
for fruitful discussions and exchanges of ideas. 

Stephane Descombes read the article in detail in the course of his first year of 
graduate studies, and spotted the defects, the typos and the errors. His criticisms 
improved considerably the article and he deserves praise and thanks for his patience. 

2. The scalar equation and the corresponding linearized operator 
We study in details the solutions of the equation 

Ut — Uyy + mu — |upM(l + ia) + = 

under the assumption a — and 

(2.1) meE\0. 
We look for solutions of the form 

If we substitute this expression into our equation, we can see that v and to satisfy 

(2.2) iLuv — Vyy + mv — + |ti|^u — 0. 

Our first and elementary results on this case are summarized in the following 
lemma: 

Lemma 2.1. Let v G Lfg^(K; C) solve ( ^.2| ) in the sense of distributions; if v does 
not vanish identically, then the following assertions hold: 

(i) V is infinitely differentiable. 

(ii) if v{y) tends to zero as \y\ tends to +oo, then lo vanishes. 
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(iii) under the assumption of (ii), the argument of v does not depend on y. 

(iv) under the assumption of (ii), the number m is strictly positive. 

(v) under the assumption of (ii), let x ~ y\frn and 

\v{y)\ ^ r{x); 

then, up to translation in space, r is the unique even positive solution of the 
ordinary differential equation 

(2.3) -mr" + mr - r'"^ + = 

which vanishes at infinity. In particular, the number m belongs to the interval 
(0,3/16). 

Proof (i) If V belongs to Lf„^(]R;C), then v" belongs to Ljo^(K;C), so that v' is 
locally absolutely continuous, and w is a function of class on R, with values in 
C. By an obvious induction argument, v is infinitely differentiable. 

(ii) Equation ( |2.2D implies that v" tends to at infinity. Therefore, v' tends to 
zero at oo: by Taylor's formula, 

v'{y) = v{y + 1) - v{y) - [ v"{y + s)(l - s) ds; 



Jo 

the right hand side of this relation tends to at infinity; therefore, the left hand 
side must also tend to 0. 

Assume w ^ 0. System (^^) is equivalent to a system of four ordinary differential 
equations of the first order in R^, and G R'' is a critical point of this system. 
We apply the theory of invariant manifolds in a neighborhood of 0; the matrix of 
the linearized system at that point has two double eigenvalues ±C, where = 
m + ioj. We make the convention that > 0. Therefore, \v\, \v'\ and \v"\ decrease 
exponentially fast at infinity. Multiply by v, and integrate over R; we obtain 
after an integration by parts 

iuj J \v\^ dy + j dx + j {m\v\^ - \v\'^ + \vf) dy = 0. 

This shows that lo vanishes, 
(iii) Define 

W{y)=m-\v{y)\^ + \v{y)\\ 
We let u = ui + iv2] then vi and V2 solve the linear differential equation 
(2.4) w" = Ww, 



and they cannot both vanish identically. If w is a solution of (2.4) which does not 



vanish identically, then w and w' cannot vanish simultaneously; if wi and W2 are two 



solutions of (2.4) which vanish at infinity, their Wronskian w'iW2 — W2W1 is constant 



and vanishes at infinity; therefore, any two solutions of (2.4) are proportional. Thus 



vi and V2 are linearly dependent, which proves the desired assertion. 



(iv) Thanks to (iii), if v is a solution of (2.2), with lu ~ 0, we may assume that 
it is real without loss of generality. Thus, it solves (2.4); moreover, v and v' cannot 
vanish simultaneously by uniqueness of solution of the Cauchy problem for ( ^.2| ). 
Therefore, we can find a continuous determination of the angle 6 such that 

(2.5) w — r cos 9,w' — r sin 9. 



It is immediate that this determination is also of class C°°. Let us multiply (2.2) 
by 2v' , and integrate, remembering that lo — and v is real; we obtain 

(2.6) + m|up - \v\*/2 + \vf/3 ^ constant. 
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Figure 2. The graph of $ for m = 11/64, to = 12/64, to = 13/64. 



By taking the hmit of the left hand side of we see that the constant on the 



right hand side of ( p.6[ ) vanishes. We use (2.5), and we can see after substituting 
into (2.6) and dividing by that 



- sin^ 



2 3 

Assume that m is strictly negative. As x tends to infinity, the upper limit of the 
left hand side of the above equation is at most equal to — min(l, |to|), while the 
right hand side of this equation tends to 0; thus, we have a contradiction. 

(v) The proof of (iv) shows that v can be taken real; define a function r by 

v{y) = zr{yy/m); 

then r is a real solution of 



—mr + mr — r + 
which decays at infinity. Moreover, if we define 

^4 







$(r) 



mr 



3 ' 



|2.6D implies 



TOT 



/1 2 



$(r). 



In order to have a non trivial solution of (2.3) which tends to at ±00, we must 
choose m so that $ can vanish for non zero values of r; letting 

^(x) = TO - ^ + ^ = x-'^Vx), 

it is immediate that "0 can vanish for positive values of X if and only if 

3 
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— (U2I^- — 



— ^JU2i^- — 



-~~ J)i)fil^'- 



"■-.-O.dMo.MB , 















— ■- — ' aii6j^_ — 








^ D-2if) 





Figure 3. The level lines of (r, s) 
m = 12/64, m = 13/64. 



^2 _ fQj. ^ ^ 11/64, 



Henceforth, we will write 



m=—{l-u) 



v = l- 



16m 



The level curves of (r, s) i—s- s'^ — $(r) are represented at Fig. ^; this figure shows 
that there exist exactly two homoclinic orbits of (2.3) through for m < 3/16; 
on one of them, r takes only positive values, and on the other one r takes only 
negative values. Therefore, we choose the homoclinic orbit which is situated in the 
half plane r > 0. The solution r is still defined only up to translation. When r 
reaches a maximum at some point xq, r'{xo) = 0, and iJj{r'^{xo)) vanishes. Therefore, 

3(1 - ^) 



(2.7) 
Since (^ 



is invariant by the reflexion x i-^ 2xq — x, we can see that r{x) = 
r{2xQ — x); if u has more than one maximum, it is periodic, which is possible only 
if u vanishes identically. Therefore, r has exactly one maximum. 

If we choose the translation parameter so that r attains its maximum at 0, then 
r is even and positive over R, which concludes the proof of the lemma. □ 

In what follows, the dependence of r or i? = over v will be emphasized from 
time to time, in which case, we will write r{x, v) or R{x, v) instead of r{x) or R{x). 
The value of R can be found explicitly; it is equal to 



1 - u 



41 



/I^ cosh 2a; 
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1/12000 

1/120000 

1/1200000 



Figure 4. The graph of R over the interval [-10, 10], for v 
lO-P/12,0 <p < 5. 



If we draw a graph of the values of r or i? for different values of ly, we can observe 
that the 'width' of r or R increases logarithmically with 1/v] see ^. 

In fact, as v decreases to i? is very close to a two front solution, and the 
distance between these two fronts can be estimated very precisely. This important 
observation in stated and proved in next Lemma: 



Lemma 2.2. Let 

(2.8) 

and 

(2.9) 



1 4 
L = - In - 

4 V 



R{x) 



V = 4e 



-4L 



4 1 



Then, we have the estimates for all x G 



(2.10) 
(2.11) 

Moreover, f 

(2.12) 

where 

(2.13) 

and 

(2.14) 



-2x 



R{x + L)~R{x) 
R'{x + L)-R'{x) < 0(l)//e-2^. 

R is a solution of the differential equations 
16 



W^--R + R', 
16 



F = -f 1 - 



4f^ 



Proof. We infer from (2^) the relation 



(2.15) 



R{x) 



3/4- R{x) 

it is also immediate that R satisfies the ordinary differential equation 
(2.16) ^' = -2^(1-4^/3), 
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which can also be found by passing to the hmit as v tends to in the differential 
equation satisfied by R: 



R' ±2Ryfi}j{R)/ m 

d a 
R{- + L 



(2.17) 

The definition (2.8) of L and a direct computation give 

(1 - i^)7?(l - 4^/3) 



(2.18) 

from which we infer the identity 

(2.19) R(- + L)-R 

1 -4i?/3 + 4z/i?V9 

Thanks to (|1|), we deduce immediately ( |2.10D from ( |2.19D . If we differenti- 
ate (2.19) with respect to x, we find 



1 - 4i?/3 + 4i/i?2/9 
i^R{l - 2i?/3)2 



R'{-+L) - R' 

R{l-2R/?,f (1 - 87?/3-f 4i?V3 



^-uR' 



(-4/3 + ^vR/%) 



(1 - 4i?/3 -t- 4t/i?2 /9)2 i_ 4^/3 _^ 4j^i^2/9 



Then, with the help of ( p. 16 ) and ( p.l5| ), relation (2.11) is clear. By a pa ss age to 



the limit as v tends to 0, or by a direct computation, we see that (|2.12 ), ( 2.13| ) 



and ( 2.14 ) hold. As x tends to +oo, R tends to and V tends to 3/16. □ 
The linearized operator around r is an unbounded operator A in L2(IR) defined 

by 

D(A) = ij2(M), Au = -mu" + Vu, V = m-Zr^ ^hr'^ . 

The spectral properties of A are important for what follows: 

Lemma 2.3. The operator A is self-adjoint; its continuous spectrum is included in 
the interval [m, -l-oo); the infimum of the spectrum of A is a strictly negative number 
X, to which corresponds an even strictly positive eigenfunction s. The function r' is 
an eigenfunction of A corresponding to the eigenvalue 0. Both A and are simple 
eigenvalues. 

Proof. Clearly, A is self-adjoint. We have 

V{±oo) = m > 0; 

The essential spectrum of yl is contained in [ni, +oo) (see Kato, [p8|, Theorem 5.26), 
and the spectrum of A is bounded from below by mina;y(a;); for all f3 > 0, the 
intersection of the spectrum of A with (— oo,m — f3] contains only eigenvalues of 
finite multiplicity; since ^ is a Schrodinger operator in one-dimensional space, these 
eigenvalues are all simple; moreover, since r is even, V is even, and the eigenfunc- 
tions of A are even or odd; in particular, an eigenfunction corresponding to the 
lowest eigenvalue A of A does not change sign, thanks to the maximum principle; 
therefore, it is even, since such an eigenfunction cannot vanish together with its 
derivative. 

Denote by s an eigenfunction of A corresponding to A; this s can be chosen so 
that 

s{x) > 0,Vx e M. 
Differentiating with respect to we obtain 

Ar' = 0. 
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It is clear that r' belongs to L^(R); therefore, it is an eigenvector of A, corresponding 
to the eigenvalue 0. On the other hand, r' is odd; therefore, we have proved that 
the infimum of the spectrum of A is a strictly negative number, A. □ 

We will obtain later much more precise information on the spectrum of A] but 
this information is easier to obtain if we take into account the nonlinearities of our 
problem ; in particular, we will asymptotically describe the eigenfunction s in terms 
of the derivative of r with respect to v, a feature which is clearly specific to our 
nonlinear situation. Thus, we have to study the dependence with respect to of a 
number of objects which appear in this section. 



(i) The Junction 

]0,1[ 



Lemma 2.4. 

u i—t r(x, v) 

is infinitely differentiable. 
(ii) The eigenvalue A(i') and the eigenprojection P{i') on Rs(i^) are infinitely dif- 
ferentiable, respectively with values in M and in the space of continuous oper- 
ators in L^CR). 



Proof, (i) Denote by H- 
belonging to iJ^(E) (resp 
icvc„(R) by 



oven (11^) 
L2 



(resp. Ll 



the subspace of even functions 



Define a mapping / from iJ, 



even( 



X (0, 1) to 



mv 



This mapping is of class C°°, we have f{r,iy) — 0, and Dif{r,v) — A, which is 
an isomorphism from -ffovcn to L^^(,jj(R). Therefore, the implicit function theorem 
applies, and r is a C°° function of v with values in H^^^^iM). 

(ii)If we divide v4 by 1 — we obtain a Schrodinger operator whose differential 
part is constant, and whose potential part depends in a C°° fashion on z/, according 
to part (i) of this proof; therefore, we can apply the results of Kato, |^ IV. 3, which 
enable us to conclude. □ 



Define three new functions (t(x, v) and a2{x, v) by 
(2.20) S{x,i^) 



aM^-AJ-^:^ = ^r{xAe-^% 
ov oL 



da{x, v) d"^ 



:r{x, 4e 



-4L\ 



The following lemma states that S + Rx , a -\-rx and CT2 + 2ax + r^x are very close to 
0; this is really a consequence of the non linearity; the benefit will be even greater 
when we will show below that a is an excellent approximation of the eigenfunction 
of A corresponding to its lowest eigenvalue. 



Lemma 2.5. 

(2.21) 

(2.22) 

and 

(2.23) 



We have the following estimates for all x € K"*", and all v £ (0; I)- 



S{x, v) + Rx{x, v) 



R{x, ly) 
dr 



< 



l-4i?(a;-L)/3 
<0(1) 



, , (9(7 , , ^ ■ / N 



< 0(1 



I' e 



THUAL-FAUVE PULSE February 1, 2008 



17 



Proof. The differentiation of ( 2.18 ) with respect to L gives 
1 



1 - 1 - 4i?/3 + 4^i?V9 



which imphes immediately (2.21). An analogous computation implies 
(2.24) A[^(^ + L,4e-4i)] = r,(a; + L, 46-*^) + (7(x + L, 46"^^), 

and 



(2.25) 92.2 



[r{x + L, 4e-4^)] = r.,^{x + L, 46"^^) 



Differentiating ( ^.18D , we find 



r(- +L,4e-*^) 



_ 2z/f(l - 4i?/3)(l - 2i?/3)2 



(l-4i?/3 + 4i^i?2/9)2 



and 



— r(. + L, 4e-4^) = - ^''^^^^ " ^^^^^^^ " " ^^/^ - 4'^^V9) 



(1 - 4i?/3 + 4j/i?2/9)3 
4i/2f(l _ 4^/3)2(1 - 2^/3)'' 
(l-47?/3 + 4z^i?2/9)4 



Observe that 

so that: 
(2.26) 



1 - 4i?/3 V 4 ^ 



,1/2 



1 - 4i?/3 
moreover we have also 
(2.27) Vx > -i 



< 0(l)^^e 



-x — L, 



1 - 4i?/3 



<0(1); 



We use the information given by ( 2.26 ) and ( 2.27 ), to infer that 



Va; > -L, 



< o(l)yI;e-^-^ 

< 0(l)y^e"--^ 



□ 



the conclusion of the lemma follows immediately. 

For later purposes, we define 

S{x + L, 4e^*^) = —^(2; + L, 4e-4^). 

From here, we will obtain an upper estimate for the lower bound of the spectrum 
of A: 



Lemma 2.6. We have 



{Act, a) 

(cr,CT) 



iv 
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Proof. Let us compute Aa by differentiating tlie equation 



-mr" + mr — r'^ + = 



with respect to L; we obtain 

(2.28) 

wfiere 

(2.29) 



Aa = p 



1-1/ 



(r^ -r'^). 



To prove our assertion, we will have to calculate (cr, a) and (ct, p). From Lemma 



2.5, we can see that 



cr'^ (x , ly) dx — 2 / rj;{x,i')'^ dx + 0{^/l'), 



and from estimate (2.11), we have 



r^{x,vfdx = 2l f'{xfdx + 0{Vv); 



but relations (2.14) and (2.16|) imply that 



Iff ( 1 



we infer from relation (2.7) that 



(2.30) 

so that 
(2.31) 

On the other hand. 



^\P'f dx = ^ + 0{V^, 



(r^ - r^ ct) = -2 / {r^ - r^)r^ dx + O(V^) 



(2.32) 



V 6 ^ 4 



Thus, we obtain from d^H) and ( |2!3^ ) the relation 



(2.33) 



(^-,a) = -^ + 0(.^/^) 



The conclusion of the lemma is a direct consequence of (2.31) and (2.33). 



□ 



Theorem 2.7. There exists a constant K such that for all small enough v the 
spectrum of A is partitioned in its intersection with the interval [—K, K] and its 
intersection with interval [2K, +oo). 



Proof. Define 



Vix) 



16 



3R + 5i?2 



It is clear that f' satisfies the Schrodinger equation 

-f " + Vf = 0. 
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As f' is strictly negative for all x, this means that is the lower bound of the 
spectrum of the operator A defined by 

D{A) ^ H^{R), Av = -w" + Vv. 

As V"(— oo) = 3/4 and V{+oo) = 3/16, the essential spectrum of A is included in 
[3/16, +oo), and thus the second eigenvalue of A is some number A > 0. 
We use now the twisting trick of E. B. Davies Q: define the operators 

~^f-d^ + Vi--L) \ + 

^ ^ -d^ + V{ L))' ^ \ A 

let ry be an infinitely differentiable function from M to [0,7r/2] which vanishes for 
X < —1 and is equal to 7r/2 for a; > 1 and let rj^ = ri{-/L). Define now a unitary 
transformation in L^(R)^ by its matrix 

jj^ _ f cos r/L - sin rjL 
^ \ sin r/L cosryL 

We compute UlJUI and we find that 



rr* _ , V'i-/L) fO -1\ d vl\-lL) fO -l\ [rf {■ / L)f 



and 

71 \ /ycos^Tyi + (3sin2 77L)/4 (cos ry^ sinTy^) (F - 3/4) \ 

\{cosf]Lsmf]L){V -3/4) V sin"^ tjl + {3 cos^ Vl) ' 

An elementary calculation shows now that 

UlJUI=J + Bl-^+Cl 
ox 

where and Cl are matrix valued functions whose norm tends to as L tends 
to infinity. It is plain that the resolvent sets of J' and A are identical. 

If E and F are Banach spaces, C{E, F) is the space of linear bounded operators 
from E to F and || the corresponding norm; ii E = F, we abbreviate 

£{E, F) to C{E). If E is the space L2(k)" ^e abbreviate C{E) to C. 

Let z belong to the resolvent set of A; then for every g in L'^{M.)'^, {J' — z)^^ 
belongs to i/2(R)2, and the function z t-^ \\{d'^ / dx'^){J — z)^^\\c(l^) is continuous 
on the resolvent set of J . 

It is equivalent to solve 

J + Bl^ + Cl 
ax 



and 



Bl^ + Cl]{J -z)-'g^ f- 



the previous considerations show that for all z in the resolvent set of we can find 
L{z) such that for all L > L{z), z is also in the resolvent set of J^. In particular, if 
we take K = A/3, we can find L such that for all L > L, the segment {K, 2K) is 
included in the resolvent set of J^. In particular, we see also that for all L > L, the 
generalized eigenspace of A relative to the interval [— min V, K] contains exactly two 
eigenvalues; the minimum of V is equal to — 21/80 + 0(i^); the above considerations 
show that these two eigenvalues tend to as tends to 0. □ 
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Let US prove now that A is very close to {Aa,<7)/{a,a). 

Denote by 11 the projection onto the sum of the eigenspaces relative to and A. 
We have 

(2.34) n = /(C - A)-' dC, 

where 7 is a circle of radius 3K/2, which is traveled once in the positive direction. 
Theorem 2.8. We have the following estimates: 

(2.35) \a-Ila\H2iR)=0{i^VL), 
and 

(2.36) A = -y +0(j.3/2). 

In particular, we can choose s = Ha as an eigenvector of A relative to the eigenvalue 
X, as soon as v is small enough. Moreover, if a is the unique solution of 

(2.37) na- = 0, Act = (1 - n) (Aa - p) 
then 

(2.38) |nCT-CT-CT|^.(j,)=0(i.2VI). 

Proof. We observe that for |C| = 3K/2, C — yl is an isomorphism from i7^(M) to 
I/^(R) which transforms even functions into even functions, and odd functions into 
odd functions. Therefore, since ct is even by construction, IIct is even too. 
We rewrite the identity Aa = p as 

P , iC-A)a 



CT = 



c c ■ 

Hence, by integration along 7, 

(2.39) na = a + ^J^{C-A)-\-'pdC. 

Since \p\ = 0{v\fV), we obtain 

(2.40) Ict-HctI = 0(iy\/I), 



If V is small enough, ( 2.4C ) implies that IIct does not vanish identically. Thus, 
by construction it is an cigcnfunction of A relative to the eigenvalue or to the 
eigenvalue A. But we know that IIct is even; thus, it has to be an eigenfunction of 
A relative to A, and 

^ _ (AHct, Hct) 
^ (nCT,nCT) 

But 

(ABct, Bct) = (p, Bct) = (p, ct) + (p, Bct - ct) = (p, ct) + 0{Lv^). 

and, on the other hand 

(Hct, Bct) = (ct, Hct) = (ct, ct) + (Hct - ct, ct) = (ct, ct) + 0{i^Vl). 

This proves ( ^.36 ). 

In order to obtain ( ^.3E| ) which holds in norm, while ( 2.35 ) holds in norm, 
we subtract the relation AHcr = AIIct from the relation Aa = p and we obtain 

"Tilf^xx — iHcr)xx \ = \p ^ AIIct — V{a — IIct)! = 0{iy\^). 
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For the last a sserti on of the theorem, we remark thanks to Cauchy's theorem, 
the last term in ( ^.39 ) can be rewritten as 

^ y^(C - A)-\-^pdC = ^ ^(C - A)-\-\i - n)pdc = a 

and a is the unique solution of 

(2.41) Ha = 0, Act = (1 - U)p. 

Subtracting {^A^) from ( |2.37D , we find that 

n(CT-CT)=0, 74(ct - ct) = A(ct - Hct). 
and this implies immediately estimate ( |2.3§| ), which concludes the proof. □ 

We need also the equation satisfied by CT2 . If we differentiate ( 2.28D with respect 
to L, we obtain 



(2.42) 
where 



Act2 + (20r3 - 6r)a^ = p^, 



dp il' , 11 s 16l^ , o, 8V , A 



dv A" ' (l-i^)2^' ' ' ' \ -v 

Therefore, 

(2.43) \p2\=0{^fLv). 
Moreover an analogous and straightforward argument shows that 

ds 



(2.44) 



dL 



(J2 



The final information we obtain are asymptotics on dX/dL on ds/dL: 

Lemma 2.9. The following estimates hold: 

dX 

dL ' 
ds 



(2.45) 
(2.46) 



dL 



6i^ + 0(i/3/2), 

= OliyVL 



Proof. For the first estimate, we differentiate the relation 

AWa = AHct 

with respect to L; we multiply scalarly the result 

dA dlia dX dUa 



by IIct and we obtain 



ina 



2anCT (dA 



dL \ dL 



— Hct, Hct 



Since 
(2.47) 

we may write 



^ = ^ {A - 3r2 + 5r4) - Bra + 20r2CT, 
dL \ — ' 



f dA \ 

( — nCT,nCTj = Av{{ir'^ -r^)a,a) + {{2Qr^ - Qr)(j'^ , a + 2a) +0{v'^L). 
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Thanks to (|2.42| ), 

{{20r^ - 6^)a^ (T + 2ct) = {Aa2 - P2, cr) 

= (ct2, p) - {P2,ct) + ((72, Act - p) + 0{v^/^); 
here we have used ( |2.37| ) and ( |2.35| ). Therefore, 



d 



—Ua.Ua] =4i.((3r2-r4)a,a) +2Aa2,(T)-— (p,a)+0(i.3/2)_ 



At this point it is clear that there exists a number c such that 
(2.48) 



but the exphcit calculation of c is tedious. It can be avoided by integrating (|2.4g 
with respect to L: we see that 

-ci//4 + 0(z^3/2^^ 



A 



and we infer from (2.36) that c = 6. 

The other asymptotics are obtained as follows: we differentiate the relation 



J 'y 



with respect to L: 



1 



dL 



2iTT 



[OA 



-dL 



We deduce from ( |2.38D that for C G 7 

\iC-A)-'a-iC-X)-' 



= o(Wl 



therefore, it suffices to find an asymptotics for 

4:V 



1 

But 



7(C-^)-^(C-A)- 



1 - 1/ 



p + {3r^ - 5r'')a + p2 - A(J2 + (C - A)a2 rfC- 



(C ~ A)-H.C - {-Aa2 + (C - A)a2) = (C - A)"! - A(C - A)-\C - A)"! 



□ 



and it is clear now that ( 2.48 ) holds. 

3. The second equation 

Let us first observe that we need only consider the case a > 0: the mapping 
a ^ —a transforms a solution into its conjugate. 

There are two different ways of writing the full problem, each with a different 
scaling. The unsealed problem is 
3 3 

(3.1) (1 - v)u" H (l-v)u-{l + m)|MpM + \u\'^u + iuu = 0. 

16 16 

Therefore, we define 

3 3 

F{u,uj,a,iy) — — — (1 — i^)u" + 7-7(1 — iy)u — (1 + iQ;)|uPw + lul^'w + iuju. 
16 16 

Whenever necessary, we will use the components of u, denoting them by ui and 
U2, so that u — ui + iu2, and we define similarly the components of F, Fi and F2\ 
therefore, we have 



F = Fi+ iF2 
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Let a = ^/s, ui — U2 — rjy^, lu = ry^. The scaled problem is 

(3.2) G(^,77,T,e,i.) =0, 
where we define 

(3.3) G{^, r], T, e, v) = Fi (^ + i^/e-q, ^fer, ^e, i^) + ^i^2 (C + i^fe r/, ^fsT, Ve, v) . 
More explicitly, 

Gi (e, r?, T, e, z.) = -rer] - ^ (1 - v)^' + ^ (1 " 

-(e-£^?)(e' + £'7') + (e' + e?7')'c, 

G2(e, r;, T, e, z.) = re - ^(1 _ ^)^" + ^(1 - z.),? 

-(e + r?)(e'+e'7') + (e'+e^')V 

Clearly, ( ^.2| ) is equivalent to (^^). We will solve the existence problem looking 
at the G formulation, and the stability problem looking at the F formulation. 
We have already determined a solution of 

Gi(e,77,r,0,z.)=0. 

it is given by ^ = r, and 77 and r are arbitrary. 

In next lemma, we state the sense in which we can find a solution for the second 
equation at £ = 0: 

Lemma 3.1. There exists a unique G M and a unique q in iJ^(R) and orthogonal 
to r which satisfy 

(3.4) G2(r,g,0,O,j^) =0 
Moreover, q is even. 

Proof. The second equation at e = can be written 

- mrj" + mr]-^^ - fr] + ^^r] = 0. 
Therefore, if we let ^ = r, we have to find q and 6 such that (pM holds, i.e. 



r — or. 



(3.5) Bq = —mq + mq — r q + r q 
Define an unbounded operator B in L^(M) by 

D{B)^H^{R), Bu = -mu" + Wu, W^m-r^+r^. 
Clearly, B is self adjoint. By construction 

Br = 0. 

The lower bound of the essential spectrum of _B is to > 0; since r is positive, 
the lower bound of the spectrum of B is zero, which is a simple isolated eigenvalue. 
Therefore, solving (3.5) will be possible if and only if its right hand side r^ — Or is 
orthogonal to keri? = M?-. This defines 6 thanks to 

J r^ dx 

Moreover, we define completely q by requiring 

(3.6) (<Z,r)-0. 

It is clear that the q we obtained belongs to i7^(R). Since the equation is 
invariant by the transformation x > ~x, q has to be even. □ 

We will need precise asymptotic information, which we proceed to give now: 
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Lemma 3.2. The number is a function of v £ (0, 1); we have the following 
asymptotics for 9 and 86 /dL: 



(3.7) 

Proof. Let us compute the difference between 3/4 and 9: 



3 

4 



/„°°(3rV4-r4)dx 



Jo r^dx 

But an explicit computation and ( 2.f 0| ) give 

3i 3, 4 



Rdx 



+ - In - + o(Vi;), 



that is 
(3.9) 

In the same fashion, 



Rdx = — + 0(^). 



/ (3i?/4 - i?2) da; = - / i?(l ~ 4i?/3) da; 
Jo 4 Jo 



7?(1 -4i^/3)da; + 0(^^) 



3 

^^(-L) + 0(V^) = ^ + 0(^^). 



Therefore, 



3 

4 



which impHes immediately (|3.7| ). 

Since r is a C°° function of € (0, 1) with values in _ff^(M), and it is integrable 
as well as cr, is a function of ly. Let us compute d9/dL; we have 



09 
dL 



(3.10) 



But 



Similarly, 



dx / 'ir^a dx - 
^0 Jo 



rfa; / 2ra dx 



r"^ dx 



Ar^adx 



16 



/ 2radx = ~ 2rr' dx + 0{^) = R{0) + 0{y^) = t + 0{y^). 
Jo Jo 4 



Therefore, the numerator of (3.10) is equal to 



— / Rdx-- R^ dx + 0(L^/V) ^ — Rn ~ -iR/i) dx + O(L^) 
16 Jo 4 Jo 16 ,' 
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128 



0{L^). 



With the help of {RM, we obtain (3 



□ 
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Let US find now asymptotics for q and related quantities. Let 
(3.11) 

r 

Then satisfies the ordinary differential equation: 

—m{(f>"r + 2<j)'r') = — Or. 

Multiplying by r we get 

(3.12) m{4>'r^)' = er^ 

Therefore, 

m{(l)'r'^){x) = a+ f {OR - R^){y) dy. 

But a must vanish: q' belongs to L'^(M) and is equal to 0r' + (fi'r; but 

r' = —rsgn{x)^/W(j)Jrn, 

hence 



(/)r' — ^qsgn{x)^/W{r)/m, 

so that (j)'r = q' — (j>r' belongs to L^(M). This implies that belongs to 
and hence a vanishes. Therefore 

1 r 1 r°° 

(3.13) mcj,'{x) = -— {9-R)Rdy = —— {e-R)Rdy. 

J-oc Jx 

Lemma 3.3. The following asymptotics hold for x > 0: 

, 31n(l - 4^(a; - L)/3) , ,- , . 

(3.14) 0'(x)--2 \ ^ J' ' +0{V^/L). 

4:LR{x - L) 

(3.15) 0(O) = |+O(1), 

o r 

(3.16) 0(i)^-_ + O(i) 

(3.17) \q\L^(^)^0{L^), Iglioedt) = 0(L), \q\^0{L^/^). 

Proof. From (|t|), we obtain 

1 /"^ 1 /"^ 

(3.18) m(/.'(x) = -— -(0-3/4) / R{y)dy~-— i3/A-R)Rdy. 

R{x) R{x) 

For X > 0, we have thanks to ( 2.10| ) 

poo poo 

/ {3/4- R)Rdy = / (3/4-^)^d2/ + 0(^^e^2^) 

^ a; ^ x — L 

= ^Rix- L) + OiV^e-^^). 



But relation ( 2.1C ) implies also that for a; > 

Rix) 



(3.19) <o(i)yi;. 

/((x — iv) 

Therefore, 

(3.20) -i- r (3/4 - R)Rdy = | + O(V^). 
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On the other hand, 
Rdy = 



i?dy + \/^e-2^0(l) 



x — L 



-\n{l + 4e^^-^''/3) +1^0(1)- 



4R{x - L) 



Therefore, we have also 
(3.21) ^ 



R{x) 



Rdy 



I- AR{x- L)/^ 
Z\n{l^AR{x-L)/i) 



R{x - L) 



+ 0(^/^). 



Thus we have obtained the asymptotic for x > 



(3.22) 



R{0 - R)dy^ R 



ln{l-AR{x-L)/3) 

R{x - L) 



Relation ( 3.14 ) is an immediate consequence of ( 3.22 ). 

In order to have an asymptotic for (^(0), we will give another asymptotic for 0', 
which will be much less precise, but more usable at this level of the discussion; 
nevertheless, the precise asymptotic ( ^.141 ) will be useful later. 

For X > L, 



ln(l -4i?(.x-L)/3) = - 



AR{x - L) 



+ R{x- LfO{\). 



Therefore, for x > L, 



<f>'ix) = -2 + ^ + 0(1):^^^^ + 0{V^). 
In the same fashion, for < x < L, 

- ln(l - AR{x - L)/3) = ln(l + e'^'-^-^^) = 2{L - x) + 0(l)e2(^-^) 



(3.23) 



<j>'{x) 



^2a;/L + 0(l/L)e2(^-''^) if < a; < L, 
^ -2 + 0(l/L) ifx>L. 

By integration in x, we get now 

^ -x^/L + 0{l), iiO<x<L, 
L-2x + 0{l) + {x-L)0{l/L), iix>L. 

The value of 0(0) is determined by the condition that q should be orthogonal to r 
(see Lemma 3.1). Therefore, 

0(0) / R{x)dx 



(3.24) <j>{x) - m 



/ (y + 0(1)^ i?(x)dx 



{-L + 2x)R{x)dx + 



Oil) 



x — L 



0{l)]R{x)dx. 



At this point we observe that 



Rdx^ — + 0{1), 



{x^/L + 0(l))i?(x) dx = — + 0{L), 



(2x- L)R{x)dx^O{L), J {0{l) + 0{l){x- L)/L)R(x)dx^O{l). 

This yields the first relation in ( |3.15 ). The second relation in ( |3.16| ) is a consequence 
of the first together with (3.24). 

Estimates (3.17) are a direct consequence of the definition ( ^.ll| ) of 0, of the 
equivalents (3.23) and of (3.15) and (3.16). □ 
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Since is an isolated eigenvalue of B for all i' in (0, 1) and since (1 — v)^^W 
is an infinitely differentiable function of v with values in the space of functions 
bounded as well as their first two derivatives in space, we can see t hat q is a C°° 
mapping from v € (0, 1) to iJ^(R), under the normalizing condition (3.6). 

Differentiation of q with respect to L will be useful too; define 

^-g(x,4e-4^). 



91 



dL 



If we differentiate (3.5) with respect to L we obtain 
(3.25) Bqx + (4r3 - 2r)qG = p + 3r^a - 9a 

where 



(3.26) 



4i/ 



P = 



{r^q - r'^q - + 9r). 



1 - V 

We recall that 9i has been defined at ( |3.8[ ) 

We will need an estimate over (}>" . From (3.12), we have 

mcj)" ^-R + 9- 2m(j>'r'/r. 

Assume that a; < and Write 

(b'{x) - 2 



then, thanks to (2.17), we have 

—2mif)'r'/r — —2(1 + a)mR' 

so that 



- R - -2m6'r'/r = 9 

4 



O(V^). 



But the following expressions are bounded independently of v for x < 0: 



'3 _ ^^1^1 _ ^ Hl-m^\-L)/3) ^^.ln(l - ARi\x\ - L)/3) . 



4 7 R{\x\ -L) R{\x\ -L) 



therefore, the term a(3/4 — i?) is an Oil/ L), and we conclude with the help of (3.7) 
that 

(3.27) \cj)"\oo=0{L-^). 

We will give now two estimates on the second lowest eigenvalue p,2 of _B; this 
eigenvalue is the lowest eigenvalue of B restricted to odd functions. 

Lemma 3.4. The second lowest eigenvalue p,2 of B satisfies 

(3.28) ^,^^ + o(L-5/2)_ 

Proof. If t2 is an eigenvalue of B, it is odd, thanks to the general Sturm-Liouville 
theory. Define an operator Bq in i^(R+) by 

D{B^)^Hl{M.+ )r\H-'{m.+ ), BoU = BuU^; 

then /i2 is the first eigenvalue of Bq. In particular, for all v G ifo(K+) \ {0}, the 
following inequality holds: 

/>oo / />oo \ —1 

H2< (m\v'f + Wv^) ( / v^dx 
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Define a test function v by 

'sin(7ra;/2(L - /L)) if < a; < L - /L, 
r{x)/r{L~ y/T). ii x > L - VL. 

For all large enough L, W < on [O, L~ VX] ; therefore an integration by parts on 
[L — VL, +oo) shows that 

(m\v \ +Wv^)dx< ^ 

^ ^ " 8(L- VI) R{L-^ 

On the other hand, 

io - 2 ' 

therefore, we see now that 

(3.29) /i2 < TT^m/Si^ + 0(L"5/2_ 

Conversely, define Li > by 



2 

then, thanks to ( |2.10| ) and ( |3.29| ) 

Li = i + iln2 + 0(l/L2); 

the function x 112 — W{x) changes sign only at a; = Li; therefore, there exists 
L2 ^ Li such that ^2 increases on [0, L2] and decreases on [L2, 00). 

Case 1: L2 < L ~ \fL. Write M — -yZ/IJ/TO; since t2 and i2 cannot vanish 
simultaneously and i2(0) vanishes, there exists a continuous determination of the 
angle 7 such that 

(3.30) t2 = 6sinM7, = 6MCOSM7, 7(0) = 0. 
A classical calculation gives 

/ ^ 2 
^-1 = ^TT^sm M7. 

There exists a constant C such that for all L large enough, 

Vx e [0, L - VI] , < -Ty(a;) < Ce-2v^; 

if we let = Cnie^^^ , 7 satisfies the differential inequality by integration of the 
inequality 

7' < 1 + £7(a:)^ 

which we integrate as: 

(3.31) 7(a;) < e"Han(ex). 

Thanks to ( |3.30| ), M"f{L2) = 7r/2; therefore thanks to ( |3.3lD , 

Af > £7r/2tan(eL2); 
but the definition of e and the condition on L2 imply immediately that 

- 2(l + Ce2)(L- VI) 
which implies immediately the desired estimate (3.28|). 
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Case 2: L2 > L — VL- We have to give a different argument: multiply the 
equation —mt'2 + Wt2 = iJ-2t2 by v and perform two integration by parts on the 
term involving a second derivative; we eventually obtain 




+ TOi2(-L- VI-0) -z/(i-\/L + 0)] 

Normalize t2 so that ^2(^2) = 1; then, by concavity on [0, L2], 

hix) >x/L2> x/Li] 
therefore, there exists C > such that 



t2Wv dx + 112 I t2V dx 



on the other hand. 



Jo 



t2Wv dx 



t2V dx > CL] 



0(Le-2^), / t2vdx = 0{VL), 



Therefore 



v'{L-y/L-0) -v'{L-Vl + 0) 



CL 



A{L-VLy 



M2 



which implies immediately ( |3.29| ) 



This proves Lemma 3.4 



□ 



4. Existence oe pulses 



We will prove that there are couples {£,,1]) which satisfy (3.2); the main idea is 
to use an ansatz suggested by previous long computations; this ansatz relies on the 
observation that for small values of e, the pulse is simply stretched; i.e., the new 
pulse obtained should be very close to a pulse with e = but a smaller value of 
v. In fact, we will work in reverse: given z/, we define — i/e~^^, with y > 0; we 
let L^ = L + y, = r(-, j/''), = q{-, v^), 6^ = 0{-, j/''), = s(-, ly^) and we have 
to determine an ansatz for e. We will start with this ansatz, with very little 
motivation. The justification will be provided by the proofs which come later. 

Once the ansatz is determined, we get an asymptotic description of the behavior 
of . Then, we go on to apply a version of the implicit function theorem cum 
estimates which is given in the appendix; for this purpose, we have to check that 
is very small, that D(^^^^^r,E)Gir\q\e\e\iy) is invertible, with 
an estimate on the norm of its inverse; this necessitates an adequate choice of 
functional spaces; finally, we check that the second derivatives of C are not too 
large in a neig hborhood of {r\q\eK,e^), and we are able to apply this implicit 
function theorem. 

Moreover, we obtain precise asymptotics on the second term of the expansion. 
These estimates are crucial for the proof of the skew stabilization. 
Define 



(4.1) 



giie)^Ci{r\q\( 
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Another expression of gi obtained by substituting the value of (— r" +r)" from ( |2.3| ) 
is 

(4.2) 51 (e) = ^^-^(r^ - r^f + e{r^q + 2r^q^ - rq^ ~ Oqf + e^{rq^ + q^f. 

1 — ly' 

The ansatz for e is given by 

(4.3) (gi(e^),s^) = 0. 

Our purpose now is to show that this ansatz determines uniquely for small 
enough values of ly, and for all y > 0. Let us introduce new notations 



(4.4) «=- ^, 

(4.5) ai^ {r^q + 2r^q^ -rq^ -eq,s), 
02 ^ {rq^ + q^,s). 

Of course, will be equal to aj with r, q, 9 and s replaced by their analogues 
q\ 0^ and 



Therefore, (4.3) is equivalent to the following equation of the second degree: 

(4.6) Kal+e^a\ + {e^fa^^ = Q. 

In order to have an asymptotic for we will need an asymptotic for aj , < j < 2. 
We introduce the notation 

0(1) 

for any quantity whose absolute value or relevant norm is bounded by some power 
of i. 

Lemma 4.1. The following asymptotics hold 

(4.7) ao = -^ + 0(V^), 

(4.8) - = ^+Ki)' 

(4.9) a2 = ^+0(L^). 
Proof. We have 

(r^ - r^ .s) = (r^ - cr) + 0{uL) = -2 / (r^ - r^y dx + 0{^/D) 







This gives (4.7). 



The second object is more difBcult. If we define 
(4.10) {r^q + 2r^q^ ~rq^ -eq,a) 



estimate (2.4C) implies that \ai — a| = iyO{l). But, thanks to (3.5) and (3.25), a 
can be rewritten as 

a = ((r2 - 9)q, a) + ^{p, q) + i(3rV -9a- d,r, q) - ^(gi, - Br). 



Thus, we have 
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We transform this expressfon using ( ^.12 ): 

(4.11) a 
and by integration by parts 
(4.12) 
But 

J ra\(j}'f dx = -2 rr' \(t)'f dx + 0{y^) 

= 1^ +2 / R4)' (t>" dx + 0{^) 

4 
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= (m(r20')'>i) - ^^("i(r20')'» ++^(/5,9), 



1 dfrt X 



\(t)' {L)f + 2 (^R-^ (t)'(j)" dx + 2 J R<j)'(t)"dx. 



From estimate (3.27) over 4>" , and the following facts 



/ |i?-3/4|da; = 0(l), / Rdx^O[l), 
Jo Jl 



we can see that 
(4.13) 



m I0T dx = 3 |(/)'(L)r /4 + 0{1/L), 
and from asymptotic ( 3.23| ), we obtain 

J ra\(l)'f dx^3 + 0{l/L). 



This impHes ( [l8[) . 

The last of the three estimates concerns 



{rq^ + q^)sdx. 

Thanks to (|3.17 ) and ( 2.35 ), the second term will contribute an 0{L^). Let us 
estimate J rq'^sdx, which will be equivalent to CL'^. More precisely, 

J rq^sdx — j rq'^a dx + i'0{\) 

POO 

= -2 rr'r'^(t>'^ dx + vO{l) 
Jo 



i?3 



i?3(0) 



4 |°° 



3 



ydx + iyO{l) 



"^^^u^ydx^ 



00 ^3 



dx + vO[l). 



From the asymptotics ( p.l4[ ) on (j)' and ( |3.15 ) and (3.16) on 0, we can see that 

/ |(i?3-i?3(O))(0'*)')| dx = 0(L3) /" |i?3-i?3(o)| rfa; = 0(L^) 
Jo Jo 

and 

/■CX3 /-OO 

/ R^\{(l)y\ dx<C a;3e-2(-~i) da; = ©(L^). 
Jl Jl 

We obtain now 



rq s dx 



1 /3\^ / 2i^ 



3 V4 



0{L 
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which is precisely the last of the three estimate of our list. 



□ 



Th eref ore, for all v small enough, and for all y > 0, the discriminant of equa- 
tion (16) is strictly positive, and there exists a unique positive solution of this 
equation, which is given by 



(4.14) 

From this relation we get 



(4.15) 



1 



^ + ^^0{\) and - = 4 + O U 

a K 4 \L 



We will apply now lemma 7.1. Define 
(4.16) C/=(e,r;,T,e), ^ {T\q\e\e''). 

Thus, we have to check that G{U^ , v) is small, and to obtain bounds for the inverse 
of D(jG{U^ , v) and for the second derivative of C? in a neighborhood of . 
The first result is a bound on G(C/^l/). 

Lemma 4.2. Assume that y G [0, L^], where p > 1 is an integer. The following 
bound holds: 

\GiU\,.)\ = KO{l). 
Proof. We have already obtained an expression for Gi{U^ , ly): it is given by 

Gi{U\iy) = K(r5 - r^f + e\r^q + 2r^q^ - rq^ - 9qf + {e^f{rq^ + q^f. 
We can see that 

1(^5 _ ^3)b| _^ |(^2g ^ 2r3q2 _ rq^ _ Qqfl + |(rq4 + q3)t| = 0(1). 



With the help of (4.15), we obtain 

(4.17) |Gl(C/^^.)| =kO(1). 
An expression for G2{U^ , v) is given by 

G2(r^9^6'^e^^/) = n{r'^q-r^ - r^q + r9f +e\2r^q^ -rq^ - q^)^ + [e^q^f 
Arguing as above, 

(4.18) \G2{U\v)\=kO{1). 



Relations ( 4.17 ) and ( 4.18 ) imply the assertion of the lemma. 

Let us compute now the differential of G with respect to J7, at (C/^, v): 



□ 



D^Gi{U\u) 

D^Gi{U\u) 
DrGi{U\v) 
D,G^{U\v) 



-m— + m + (5r4 - Sr^)^ + \e{2rq - q^ + 6r^q^) + e^q^] 
dx^ 

i^A^ + «(5^4 _ 3^2)b _^ [^(2rg -q^ + 6r2q2) + s^q^]\ 
[e{r^ ~2rq-0 + Ar^q) + e^{3q^ + Arq^)]\ 

-[£q]\ 

[-eq + r^q - rq^ + 2r^q^ + e{2q^ + 2rq^)]\ 
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and, similarly, 

D^G2{U\ v)=[e- - 2rq + Ar^q + e{-q^ + Arq^)\\ 

D^G2{U\ v) = ^—^B' + K{r^ - r^f + [e{6r^q^ - 3q^ ~ 2rq) + ?>e\'']\ 

DrG2{U\iy)=r\ 

D,G2{U\ I/) = [-rq^ -q^ + 2r^q^ + 2eq^]\ 
Therefore, we may write 

DuG{U\u) = {I - K)g'' + kH, 

where 



w B r z 

and the letters w, z and z denote multiplication operators by the functions 

w = e - ir^ - 2rq + ir^q, 

z = -9q + r'^q - rq^ + 2r'^q^ 

z = -rq^ -q^ + 2r^q^. 
The matrix H has components Hij, 1 < i < 2,1 < j,< 4 which satisfy 

max|i7yU = 

Lemma 4.3. Let E2 be the space 

E2 = {^e Hl^^iR),^ ±s}x {77 e Hl^,M,V ± r} X 
and let Eq be the space 

Eo = 

The space E\ is defined in an obvious way. Then, for all p < 00, for all small 
enough v and for y < , D(jG{U^ is an isomorphism from E\ to Eq, and the 
norm of its inverse is bounded by 0(1). 

Proof. We first study Q as an operator from E2 to Eq. Given h = {hi, /i2)"^ G Eq, 
we first solve 

(4.19) gui ^ h. 



(4.20) 



where Ui = (^1, 771, ti, ei)'^ is to be found in E2. Equation (4.18) is equivalent to 
the system 

'A^i + zei = hi, 
^w^i + Brji + rri + z£i = /i2. 
Since ^1 is perpendicular to s, the first equation gives 
(4.21) {z,s)ei^[hi,s). 

But {z, s) is p recis ely equal to ai defined at ( [4.5| ). Therefore, there exists a unique 
El satisfying ( 4.21 ) and 

\ei\<0{l)\h\. 

Restricted to even functions perpendicular to s, A is invertible, and its inverse is 
bounded thanks to Theorem 2.7. Therefore 

161 = 0(1), 
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and by classical elliptic estimates 



For the second equation of (4.19) to have a solution, the following orthogonality 
condition must be satisfied: 

/i2 — zei — w^i — rri _L r. 

Hence 

_ (/i-2 - zei - w$i,r) 

|,|2 

which implies immediately that 

Inl < WHeM^). 

The operator B restricted to r-^ is invertible, and the norm of its inverse is estimated 
by = 0{L^), thanks to Lemma 3.4. Therefore, the unique r]i orthogonal to r 
which solves the second equation of ( 4.19| ) satisfies 

\m\ < IMEoOii). 

By classical elliptic estimates 

\m\HHM) < \\h\\EoO{l), 

and therefore 
The full system 
can be rewritten 



||C/l|b. < \\h\\EoO{l). 

(1 - K)g^Ui + kUUi = h 



1 — K 1 — K 



Clearly, 



\\{Q'r'nciE^,)^o{i). 

Therefore, for v small enough and y < L^, we can see that 



1 - AC 

and our assertion is proved. □ 

We have to estimate D^G{U, e, v) in a neighborhood of (C/^, . Let us choose 
for this neighborhood the ball of radius 1 around [/'', in the space £"2- Without any 
calculation, we can see that D^G contains only multiplications by functions which 
are polynomials in ^, 77, r, e and v. Therefore, the norm of this operator is 0(1) 
for y < LP. 

Now, these results enable us to apply Lemma 7.1. There exists for all small 
enough ly a unique U £ E2 such that 

G{U, v) = 0. 

Moreover, 

■\\U-U'\\e^^=kO{1), 

\\U-U^ + DuG{U\v)~^G(Jj\ = K^O{l). 

It will be useful in what follows to have an asymptotic for 
Ui ^ DuG{U\v)-^G{U\v). 



(4.22) 
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The information needed is summarized in Lemma 4.4. But we need to introduce 
new notations. The orthogonal projections P on s, Pi on s^, Q on r and Q± on 
r'^ are defined by 

Pit; = u — Pt;, Qi_v = V — Qv 

We will also abuse notations: when / is orthogonal to r' we denote by A~^f the 
unique solution u of Au = f which is orthogonal to r'; the analogous convention 
will be used for B. 

Lemma 4.4. Let 

The following asymptotics hold 

(4.23) ^i = {A^)-^Gi{U\i^) + n^O{l) 

(4.24) 771 = {B')-^Q'^{G2{U\iy) - w'^i) + k^O{1), 



(4.25) 
(4.26) 



(G2(^/^^^)-$l7«^r^) 



rb\2 



ei = K''0(1). 
Proof. The vector f/i satisfies the equation 

(4.27) g^u_ 



1 - K 



1 - K 



The first equation of (4.27) yields 

1 — K 

Since HC/iH^b — kO{1), we can see that ( 4.26| ) holds. Therefore, 

A^^i = Gi{U\i^) + K^O{l), 



which together with the orthogonality i^i _L implies (4.23). In a similar fashion, 

G2{U\iy) K 



+P'ryi +rVi 



1 - K 



1 - K 



(Ht/l)2, 



so that 



G2{U\^) 
I- K 



Hence ( 4.25 ) holds, and ( 4.24 ) is inferred immediately from ( 4.25 ). 



□ 



The only remaining task is to check that the (^, 77) obtained in this fashion yields 
an u = ^ + i^T] which has the adequate asymptotic behavior at 2; = cxd. The 
linearized equation at infinity for u is 

3 3 

(1 - v)u" H (1 - v)u + i^/eru = 0. 

16 16 

Therefore, if [3 is the square root of (l + iy/er'^ /m which has positive real part, 
there is a constant 7 such that 

The phase of u is given by 



argu(x) — arctan- 
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According to our ansatz, this is very close to 

arctan(^/£(/)''(x)), 

whose behavior is exactly the behavior we postulated at the beginning of this article. 
Finally, we obtain: 

Theorem 4.5. For all p > 0, for all v small enough, there exists C > such that 
for all a satisfying 

(4.28) \a\<a„, = ^V^{l-CL-P), 

there exists a pulse solution of ( [l.7| ), i.e. a solution of the form 

u{x,t) = e'"^V(a;)e*'^(^\ 

where r is a positive function which decays exponentially at infinity, and 4>{x) is 
asymptotic at infinity to —C\x\ + D, with C a positive constant and D a real con- 
stant. 

Proof. We just have to translate the condition on a from the condition on We 



know from (4.15) that there exist C", k and K such that 



- \4 L + yJ 

Since the function 

1 _ c" \ 1 - e-'^y 

2 L + y) I- Aue-'^y 
is increasing on IR+ for v small enough, we can see that if 

,2 ^ M ^ \ ., ^ ^ .,2 T^Tk 



y ' 



(4.29) < [ ly n- - KL", 

the existence theory works; but it is clear that for i/ small enough, the right hand 



side of (4.29) is at least equal to z^(l — C/L^)^/4, and the theorem is proved. □ 



5. Stability of the pulse 
Our purpose is to study now the spectrum of 

DuF{£, + iVsiJ, t, y/e, v), 

under the assumption 

the parameter v is small enough and there exists p>l such that < y < 
We define 

(5.1) V = DuF{^ + i^er^,^reT, ^fe, ~ k), 
and we observe that 

(5.2) V^A^+^B + kC, 
where 

(...) A^^ii IX s-d ''■^)»ac.('^i^ I 
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More precisely, we have 



(5.4) 
(5.5) 
(5.6) 
(5.7) 



Sl2 
621 
Cll 
C22 



V^(t - - 2^77 + 4^3/;) + £3/2 (-r/2 



4er;3) 



5^^ 



^(1 - k) 
+e(2e?7-r;2 + 6^^772) 



1 - K 



(5r4 - 3r^f 



1 



(r^ 



1 
K 
1 



It is immediate from ( [4.22D , ( [4.23D , (|^), ( |4.25| ) and ( |4.26| ) that the foUowing 
estimate holds: 



12\\L° 



IS21I 



22\\L=°( 



kO{1). 



Identifying the multiplication by S or C to an operator from L'^{K.)'^ to itself, we 
see that 

(5.8) \\B\\c + \\C\\c = 0{l). 

The first result concerns the dimension of the generalized eigenspace associated 
to the eigenvalues of 2? contained inside a conveniently small circle; for this purpose, 
we let Ki > be a constant such that for all L large enough 

Ai2 + A > 4/{KiL^) 



Theorem 2.8 and Lemma |3.4| guarantee the existence of such a Ki. 

Lemma 5.1. Assume (|^); let 7 be the circle of radius 2/ {KiL^p^ and of center 0, 
traveled once in the positive direction. Then, the spectrum of T) does not intersect 
7; the eigenprojection associated to the part of the spectrum contained inside 7 is 
of rank 3. 



Proof. From Lemma 2.6 and Theorem 2.8, we know that the spectrum of A inside 
7 contains only the semisimple double eigenvalue and the simple eigenvalue A*". 
We know also that is the closest element of the remainder of the spectrum of 
to 7. Define 

d(C)=niin(|C|,|C-A^MC-M2|)- 
We define a twice punctured disk D by 

i?^ = {CeC: Id <M^,Cy^O,Cy^A^}. 

Then we have the following estimate on the resolvent TZ'^{C,) = (C ^ -^'')^^ oi Ai: 

vce^?^ ||7^0(c)||^<lMC). 

It will be convenient to embed T) into a holomorphic family of operators "Die) 
defined by 

V{c) =A^ + cB^- c'^C, 

and to do the same type of calculations as in Chapter II of ||2^ ; however, the results 
quoted therein cannot be directly applied, since our eigenvalues are not uniformly 
isolated with respect to the remainder of the spectrum. 

The Neumann perturbation series for ^{c) is the expression 



(5.9) 
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this series converges provided that 
(5.10) 



c| <min(d(C)/(||S||^ + ||C||^),l) 



We infer from (5^) that there exist K2 > and ^2 such that 



\\B\\ 



\C\\r <K.L''': 



without loss of generality, we may assume that for all the L's that we consider, 
K2L^'^ is larger than 1; thus the function TZ is holomorphic with values in C on the 
set 

(5.11) A = {(C,c) :CeI?^L'=^if2|c| <d(C)}, 
and the following estimate holds over A: 

(5.12) \\n{C,c)\\^< 



1 



d(C) - |c| K2L'^- ■ 

If 7 is the circle defined in the statement of the lemma, any element of 7 satisfies 
for L large enough d{C) > 1/{KiL'^p). Therefore, the circle 7 is included in the 
resolvent set of 2?(c) provided that 

KiK2L^P+''^ |c| < 1. 

This will hold for c — k and L large enough, and the first assertion of the lemma is 
proved. 

Since TZ is holomorphic in a neighborhood of c = 0, it admits a Taylor series of 
the form 



7^(C,c) = ^7^^(C)c'", 



J=0 



with 



7^^(C) = ^ J nic,c)c-^-'dc, 



where 7' is a circle about of radius strictly less than d{() / (^K2L''^) . The first 
term of the expansion is equal to 

(5.13) 7^0(C)-(C--4^)"^ 

hence our notations are coherent. If we take 7' of radius d{()/ (2/^2^*^^), we have 
the estimates 



(5.14) 



- d{c) V d{0 

and also, for all c such that |c| < d{(^)/(2K2L'^^), and for all fc > 0, 

1 



(5.15) 



fc-i 

3=0 



< 



2i^2i'^icr' 



d{C) - 2K2L''^\c\ 



The projection on the total eigenspace of T>{c) relative to the eigenvalues inside 
7 is given by 



-Pic)^^ f 7^(c,c)dC• 

2«7r J 



We infer from ( 5.14 ) that the operator defined by 

1 



2iTT 



n^OdC 
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satisfies the estimate 



for all c such that 2KiK2L^^+'^''\c\ < 1 and for all fc > 0, we have 



(5.16) 



V{c) 



k-l 

3=0 



< 



A{2KiK2L^^+^P\c\f 
1 - 2^KiK2L''^+^P ' 



For 2 KiK 2L''^'^'^p\c\ < 1, 7^(c) is a holomorphic function of c. For c = 0, rela- 
tion ( ^.13 ) implies that Vq is the projection on the space spanned by s, ir and r'; 
therefore, it is of dimension 3. A classical argument shows that a continuous family 
of projections which is of finite rank for some value of the parameter is of finite 
rank for all its values. This proves the second assertion of the lemma. □ 

Let us define now a basis of Ini7'(c) which depends analytically on c in a neigh- 
borhood of 0; since we expect to find a nontrivial Jordan block, we will not look 
for a basis of eigenvectors of I?(c). 



Lemma 5.2. Define vectors 



(5.17) 



w{c) = V{c) 



v{c) = V{c) 



t{c) = V{c) 



(r^)' + c2(e'- (r^)')/K 



There exist positive numbers > 2K1K2, > k2 + 2p, and such that for 
K^L'^^IcI < 1 the vectors w{c), t{c) and v{c) are holomorphic in c and form a basis 
oflmV{c); in this basis, the restriction of the operator V^c) to Im7^(c) has matrix 
M{c); M{c) is analytic in c and for KjL'^^lcl < 1, it admits the expansion 

00 

M(c) =^MV, 
3=0 

with the estimate 

(5.18) ||^'IL(C3) <i^4i''^(2i^3i'-')'- 

Similarly, w, t and v admit expansions 

00 00 00 

w(c) = cP , t{c) = t-' c' , v(c) = v^c-' 

3=0 j=0 3=0 



with the estimates 



max \w' 



Finally, when c ^ 



< KiL''^ {2K3L''^y . 

^K, M(c) has the form 

'Afii(v^) o^ 

M{^) = ( Af2i(V^) 
0; 

Proof. It is clear that w, t and v are holomorphic functions of c if 2KiK2L'^'^^'^p\c\ < 
1. By construction, they belong to Im'P(c). For c = 0, w(0), t{0) and v{0) are 
orthogonal vectors spanning 7^(0); therefore, there is a neighborhood of c = for 
which w{c), t{c) and v{c) constitute a basis of Im7^(c). Let (, ) denote the C-bilinear 
product on L^(R;C) defined as follows: if 



X 
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belong to L'^{R; C)^, we let 

{X,Y) = / {Xi{x)Yi{x) + X2ix)Y2{x)) dx. 



In order to estimate the size of the neighborhood of where w{c), t(c) and v{c) 
constitute a basis of lmV{c), we introduce a quasi Gram matrix r(c) defined by 

/{w{c),w{c)) {w{c),t{c)) 
r(c)= {tic),wic)) (t(c),t(c)) 

V {vic),v{c))^ 

The terms {w{c),v{c)) and (t(c),u(c)) vanish because w and t are even functions, 
while V is odd. We also remark that r(c) is symmetric, not hermitian, and that 

/3/8 + 0(\/I>) 

r(0) = 3(i + 2;)/2 + 0(l) 

V 3/8 + 0(\/j>), 

We infer from the hypothesis |c| < 1 and from (4.22) that there exist > 2K1K2 
and > 2^2 + p such that 

(5.19) ||r(o)-i(r(c)-r(o)||^(„3) < \c\k,l'''. 

Therefore, if |c| K^L*'^ < 1, r(c) is invertible, its inverse is analytical in c and it 
satisfies the estimate 

|r(o)-M 



|r(c)-i,,^3, < 



l£(C3) - i_ |c|if3Lfc3 ■ 
Denote by 

He) tic) vie)] 

the function from R to 3 x 3 complex matrices whose columns vectors are respec- 
tively equal to u'(c), t(c) and w(c). We define three vectors wic), i(c) and vie) 

by 

[wie) tie) vie)] = [Vie)wic) 7'(c)i(c) ■P(c)w(c)] r(c)^i 

Then, for K^L'^^lc] < 1, liic), tie) and u(c) are holomorphic in c and they constitute 
a basis of ImP(c)^ which is dual to the basis w(c), tie), vie) of Im7-'(c), i.e. 

[wic) tic) w(c)]"^ [?I>(c) tic) vie)] dx = lc3. 



Thus, Af(c) is given by 

Mic)= j [wie) tie) {)(c)]^ [2?(c)u;(c) 27(c)i(c) I?(c)v(c)] dx. 

It is immediate that I?(c)ti;(c), I?(c)t(c) and I?(c)?;(c) are analytical in c with values 
in L^(R; C)^: in fact, they are polynomial in c. The study of section ^ and in par- 
ticular estimates ( f.22| ), completed by estimate ( |5.1(]| ) for fc = and estimate ( 5.19| ) 
show that there exist numbers Kj^ and fc4 such that the following estimates hold: 

if |c| = l/(2i^3L''^), then 

max(|w;(c)| , |t(c)| , \vie)\ , \wie)\ , |£(c)| , |{)(c)|) < Kih'^- 

and max \M,kic)\ < K^L'"'^ . 

l<i,j<3 



Standard methods of analytic functions give (5.18) and analogous estimates for y, 
t and V. Finally, differentiating with respect to (3 the relations 

Fie'^{^ + i^i^),^r,^,v)=QFii£,+i^T^i-+l3),^eT,^e,v)=Q 
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we find that and are eigenvectors of I?(v^) relatively to the eigen- 

value 0. Thus w and t belong to ImV. Therefore the form of the matrix M(y^) 
is clear, and this concludes the proof of the lemma. □ 

All our effort will be concentrated no w on getting the first three terms of the 



expansion of A/ii(c) since we infer from ( 5.18 ) that 



|Mii(Vk) - Mn - VkMI^ ~ kMIA < k^/^0{1). 



Lemma 5.3. The first terms of the expansion of V have the form 
(5.20) P.= f° ^k). r--'"" » 



and the first terms of the expansion of w and t are given by 

»'=(:)■ '"=(»)■ '"^C° 

with 

(5.22) si = (A^ - B^)-^Q^^B2is^ and n = -{A^)-^ P^^B^r^ - ^I^P^. 

V K 

Proof. The first statement of the lemma is an immediate consequence of the defi- 
nition of the V^^s and of the special form of , B and C. Let us calculate 7^2i^^- 
7^l(C) is equal to 7^0(C)67^°(C) and thus 







,0; " \{\' ~0-\C-B')-^B2is' 
Hence, thanks to the theorem of residues, 

si = Vks' ^^jiC- B'r'B2is'^^ = (A^ - B'r'Q'^B^is'. 
Similarly, according to the definition of t{c), 



and 



Vl^r' = 1 / (C- A^)-i6i2r^^ = i^A'r'p'^Bur'. 



□ 



Now we are able to give expressions for the first coefficients of the power expan- 
sion of Mii(c): 

Lemma 5.4. The coefficients M^^, for < j < 2 and Mi^-^ for j = 0,1 are given 
by 

(5.23) Mi^i = X\ M^, = 0, 

(^9A\ n A/1 iB2is' -X's,,r') 

(5.24) Mil = 0> ^^21 = pp > 

,.2 iBi2Si+Cns' -M^,r,,s') 

(5.25) Mil = TTT^ ■ 
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Proof. From the analyticity properties we infer that 
(5.26) ^^° = AfiV + M2"it°, 



(5.27) 



(5.28) 



2 



From ( 5.26| ) and the value of t{0) we infer immediately ( ^.23| ). Relation (5.24) can 
be rewritten with the help of (5.21) as 











Ml, 



Ml, 



It is immediate that Ml, vanishes. Moreover, we perform the scalar product of the 
second component of the above identity with ; since 
we obtain t he se cond relation of ( 5.24 ). 
Finally, ( 5.28 ) implies that 

A^S2 + Si2Si + Ciis^ = A^S2 + Mns^ + Ml,n. 

If we perform the scalar product of this relation with s' and if we observe that 
{A^S2,s^) = (S2, A^s^), we obtain ( |5^ ). □ 

We give now an asymptotic of Ml, in terms of r^, and a^: 
Theorem 5.5. The coefficient Mf, has the following asymptotics 



(5.29) 



K 



da 



dL 



where a has been defined at (4.1C). 

Proof. According to ( p. 35 ) and Lemma 4.4, we have 

^215^ = ^/^[{9 - 3r2 - 2rq + 4r^q)a] ^ + (k + i^^)O(l), 
and in virtue of ( 3.25| ) and ( [3.2(]| ), we can see that 



(5.30) 
Therefore 
(5.31) 
and 



^215^ 



Si 



K 



-qU + (« + ^')o(i) 



Ml,^-^e\ + iK + iy')o{i). 



Similarly, according to Lemma 4.4 

^.2.^ = /fM 



Bi 



+ r'' - 2rq + 4:r-\)r\ + kO{1), 



and thanks to (3.5) we have 
Bur^ 
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Therefore 



and in particular, 
(5.32) 

There remains to calculate 



{B12S1 + discs'' 



We infer from (5.31) that 



Bi2\j-Q'qls' ) +(k. + j/^)0(1) 



and we have the following relation 

iB,2r\s') = .[^{{Aqf + nOil), s') = (k + ,y')0{l). 

Therefore 



(S12S1,/) = \l-iBi2qls'') + iK + i,')0{l) 



-{{-e + - 2rq + 4r^q)qi,a) + {k + i^^)0{l). 



On the other hand. 



{Cns\s') = ((5r4 - 3r')a, a)' + -{[5^ - 3^^ - {5r' ~ 3r'f]a\a') 
+ -{{2rq - q^ + 6r^q^)(T,af + kO{1). 

Hi 

But according to ( 4.22D , a Taylor expansion gives 

5^4 - 3^2 - (5r4 - 3r2)^ = {20r^ - 6rf{^ - r^) + ^^0(1) 



and thanks t o Lem ma 4.4 , this expression is equal to —{20r^ — 6r)^i + k^O{1). 
Hence, from (|2.42|) 

(5.33) 

= (4,Gl(C/^^.))+«^.^0(l) 
= (ct^, K(r5 - r^f + e{r^q + 2r^q ~ rq^ - Oq^) 
+ k{v^ + k)0{1). 

Now we obtain the expansion 
(Cns^s^) + (Sl2Sl,s^) 

= ((5r4 - 3r2)a, a)^ + -{{-6 + - 2rq + 4r3g)qi, a)^ + (r^ - r^, (73)^ 

+ -(r^q + 2r3q2 _ _ ^^^b ^ £^(2r<z - q^ + 6r^q^)a, af + {fi + J^^)O(l). 
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We observe that 

((5r4-3r2)CT,a) +(a2,r5-r3) 



{r^ _ r^)(r^.^^ + 2a^) dx + 2 rl{5r^ - Sr^) dx + V>^0{1) 



therefore, with the help of ( 5.32| ), we are left with 



d 

— ir^q + 2r'^q^~rq''-eq,a) 



This concludes the proof of the theorem. □ 

Therefore, the only remaining question is to obtain an asymptotic for dai/dL. 

Theorem 5.6. The following asymptotics holds: 

da 37r2 



dL 32L2 



Proof. We recall (4.12); we expect that the dominant term in a will be 

d 

oL 

the remaining terms being small with respect to b, as we will check later. The 
number b is also equal to 



b ~ m J ((T^ + r(T2) |0'| dx + 2m J racj)' cjy'i dx. 
Thanks to estimate (2.23), 

J (^'+ra2)|0f dx^2j^^[a^+r{-2a'-r")] \cj)'f dx + 0{V^), 



and by integration by parts: 



{a^ + ra2) \<j>'\' dx ^ 2 I [{cr + r'f \(f)'\' + 2r{2cr + r')<j)'(t)'^igr] dx + 0(^^) 



Now, thanks to (2.22), we obtain 

b = 4m ra4)'{(t)" + dx + ^7^70(1). 



We differentiate ( 3.13| ) with respect to L, and we recall the definition ( 2.20 ) of 
S\ then 



,,dm 1 
b 1 

dL m 



ft 

{S + R') {e-R)Rdy-—— {9 - 2R){S + R') dx 
Jx K[x) 



In the above expression, the principal term is the first term on the right hand side; 
the other terms are estimated as follows: it is immediate that 

(t)'dm/dL^ I'Oil); 

thanks to ( ^.2lD and ( |3.20D , 

j {e-R)Rdx^O{V^); 
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thanks to ( ^.2lD , ( |3.19| ) and ( |3.2l| ), we also have 



therefore 



^ r {e-2R){S~^R')dx < / 



m 



- L) J.-L 1 - 4i?(y)/3 

/•OO 

iy Rdy + 0{V^). 



dy - 0(i.); 



i? 



Thanks to and ( ]2^ ) 



ar = -^R'{- -L) + 0(l)i.3/4e-2(-i); 



it is straightforward that 



Jo R(x^L) ^ 



and 



/>oo 

/ 0(v^)(/)Vcrda: = 0(v^). 
Jo 

Thus, we have proved that 

r ^^ln(l-4^(•-L)/3) 



R'{-~L)- 



R{- - L) 



dx + 0{Lv^/'^). 



In order to give an asymptotic of the integral in the above expression, we cut it 
into three pieces: one piece from to L — \/Z, which is 



the second piece is 



l-Vl R{—L) 

which we integrate thanks to the change of variable y = AR{x — L)/3: it is thus 
equal to 

JO y 

The last piece is 

L-VL R{- - L) 

since for x > L — \/L, <f)' + 2 ~ 0(l/vX), this last piece is an 0(1//L). Finally, 
we have obtained 

There remains to estimate the other terms in da/dL. It is easy to see that all of 
them are of order i/0{l), and therefore negligible before the error given in the above 
formula. This completes the proof of the Theorem. □ 



We are now in grade to state the stabilization property: 
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Theorem 5.7. For < y < , the following asymptotic holds: 



+ KO{L-'''^) + {v^f'^0{l). 



(5.34) Mn(V^)=-^-, 

There is a critical yc such that aijjc, v) vanishes; an equivalent for is given by 

1 



Vc ^ 

The corresponding critical parameter is 
(5.35) ' ^ ^ 



InL. 



0(L-5/2) 



2 ' ' V' 

Proof. The asymptotic ( ^.34 ) is a consequence of all the previous results. In order 
to find the critical value of y, we have to solve the equation 



en 



2 ^ A{L + yY 



After replacing n and by their respective values, and defining the new unknown 

Y = e"4^, 

we find the following equation in Y: 

n^l-Y) 



Y 



O{L-'/')=0 



24(L- lnr/4)2 

It is clear that there is a solution Yc of this equation which satisfies 



Yr, = 



2AL' 



+ 0(L-5/2). 



Using relations (gj), ( p^ , ( p2| ) and (|3|), we can see now that at the critical 
value of 2/, 

^2 



1 - 
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Relation (5.35) is an immediate consequence of the above, together with a 



The final result of this section is 



□ 



Proposition 5.8. Let u be the solution defined at theorem 4-5. Let \a\ < am be 

the number defined at (4.28) and ac be the number defined at (5.35). Then, the 
solution u is stable iff en c < \ot\ < am, and unstable if \a\ < Uc. More precisely, 
in the first case, the spectrum of the linearized operator at u contains exactly one 
negative eigenvalue, if \a\ < ac; when alpha — ac, the eigenvalue is of alge- 
braic multiplicity 3 and geometric multiplicity 2, with a non trivial Jordan block of 
dimension 2; in the second case, the linearized operator at u has the semisimple 
double eigenvalue 0, and the remainder of the spectrum is included in Sftz > /i > 0. 

Proof. This Proposition summarizes our previous analysis: the part of the spectrum 
of T> which is outside of the disk of radius L^'^'p^^ around is contained in the half 
plane 

SRC > CL-^P, 

Lemma 'iA and the definition ( |5.2| ) 



for V small enough, according to Theorem 2. 



of V. Therefore, the change of stability is equivalent to the change of sign of 
Mii(y^) given by ( ^.34 ). Theorem 5.7 implies that Mn vanishes for some value 
Uc of |a|, whose asymptotic is given by ( 5.35| ). The results on the multiplicity of 
will be clear provided that we show that A/21 does not vanish; but formula (5.3C) 
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together with (3^) proves that this is the case. This concludes the proof of our 



proposition. □ 

The last step is to understand the evolution of small perturbations of the pulse; 
the linearized evolution is easy to understand: if 6i{t) are the coordinates in the 
basis {z,iu,u'} of the projection of the linearized perturbation on the space ImQ, 
then 

g-at 





For a > 0, we can see that the asymptotic behavior of the perturbation is simply 
described by a change of phase. The smaller a, the larger this change of phase, and 
the smaller the convergence to the phase shift. 

This linear theory can be justified, using Chapter 5, Exercise 6. Let us 
denote by u {x, t ) the stable pulse that we have found, and by u{x, t) the perturbed 



solution of (1/7). It is clear that S is a sectorial operator, because it is a bounded 
perturbation of a positive self-adjoint operator. Then there exists in the stable case 
a i5 > such that if \uo — u{-, 0)\h^(u) < then, there exist X{uo) and vE'(wo) such 
that 

lim |u(a;,i)-u(a;-X(uo),i)e^*("")|ff2(R) ^0. 

Moreover, a straightforward calculation shows that (5) gives the correct behavior 
for small S. 

6. The general situation 
It is indeed possible to perform the same kind of computations for the general 



case, i.e. equation (1.8), under the assumption 

(6.1) '^m.j = ^/efij, j = 1,2,3. 

let us give the main steps of the calculation; the functional validation is entirely 
identical, and we will shorten the exposition by giving only the relevant part of the 
expansion, and without bothering to give estimates on the remainders. 



Under assumption (6.1), the functional F becomes 



|4 

u. 



F{u) = iuju — mu" + (m + iy/e iii)u — (l + j\/e/^2)|upM + + ^V^A^s) 1^1 

The functions Gi and G2, which are analogous to the ones defined at the begin- 
ning of section ^ are given by 

Gi it V, T, £, 1^) = ~£TT] - m^" + - fuer] - (^^ + ejy^) {£_ - H2£v) 

The second equation ( |3.5| ) is changed into 

Bq = -Or - nir + ^i2r-^ - tJ^^r^ , {q,r) = 0, 
and 6 is determined by the orthogonality condition 

J ((0 + Mi)i?-M2i?' + Ai3i?')- 
For (f) given by ( 3.11 ) , we have the relations 
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and 



m(j)'{x) 



1 f°° 

— - / [-OR - fiiR + /i2i?2 - /^3i?^] dy. 



The derivative qi of q with respect to L satisfies the relation 

Bqi = -9ir -6(7 - fiia + 3^27"^ cr - 5^i3r'^a + 2rq(j - Ar^qa + 
The asymptotics for 9, 9i and 0' are given by 



la 
8i2 



3^2 27/^3 
8L 64L 



y^o(i), 



3^2 _ 27/i3 
64^2 



and 

(6.2) 



^ ^(1 + e2(--^)) ln(l + e2(^--)) + C2 + ^(l 



=2(2;-L)^ 



where the numbers ci, C2 and C3 are given by 
(6.3) 



3At3 



9^3 „ , 3^*3 

Cl=/i2 ^, C2 = -2fl2 + —, ^3 = 



We obtain these asymptotics with the help of the following calculations: 



R{y) dy = -ln{l + e 
Riy^dy ^^[\n{l + 



-2x\ 



-2x\ 



„2£C\ 



/ 

-/ X 

f dy = ^^ [ln(l + e-'^) - (l + e^'^' - (l + e'^" /^] ■ 

The existence ansatz is completely analogous to the one defined at the beginning 
of section ^ by ( [4.1[ ), ( p^ ) and (4.3); the definition of aq is as in ( [4.7D , and we let 

ai = ((-6* - Ail + M2r^ - ^l:^r% - q^r + 2r^q^, s) 
a2 = (^29'^ + rq'^ - 2^3r^g^, s), 
03 = -(m39^,s)- 

It is straightforward to prove that the new ai is still equivalent to ■m(r(T(p' , <j>'^ 
as in ( i.l2 ), and we jus t ha ve to calculate (p'iL) in order to apply formula ( fl.l3| ); 
we infer from estimate ( |6.2| ) that 



Therefore, 
(6.4) 



ai - ^(-2Ai2 + 15Ai3/8)' + 0{1/L). 



:(-2^2 + IS^is 



The remainder of the existence proof is identical to the proof given in section 
its details are left to the reader. 

Next step is to calculate DuF at the solution obtained by this existence proof; 
we will have (5.2), with B and C of the same form as in (5.3), and the new values 
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for the coefficients of cCilb and C are given by 



6 



12 



S21 - 



Vk{1 - k) 



C22 — 



V^(l - k) 

1 r -|t/p - 2e + 2£^2e^ + 1^1^ + 4|Mpg^ - 4A^3e^|M|'£ 

Kl 1 — K 

1 r-|Mp - 2r]'^£-2r]^fi2e+ \u\^ + 4\u\'^ri'^e + ifj.3\u\^T]^ 

Kl 1 — K 

The principal part of B12 is 

^(-6' - Ml - 2rg + /i2r2 + Ar^q - fisr^f, 
and the principal part of B21 is 

.[^{e + Ml - 3^2^' - 2rq + Ar'^q + b^isr^f . 
The principal part of Cn is 



+ (3r2 - 5r4)^ 



('2 4\ ^ 

r — r 



Lemmas |5.l|, |5.2| and 5.3 have an identical statement and an identical proof. The 



calculations of Lemma |5.4| are modified only to get Mfi : we have first 



+ 



With the same computation as in (5.33), the last term of the principal part is 
((72, Gi{U^ , i/)) /k which is equal to 



((72, - 



(72, -6q — fiiq — rq + ^27" q + 2r q — fi3r q 



up to higher order terms. On the other hand 



(Si2Si,s^) ~ .[^{Bi2q\,s^) ^ -{i~e ~ fii - 2rq + fi2r^ + ir'q - ^i3r'')qi,ay 



We find finally the same result as in (5.29) 



The computation of da/dL proceeds along the lines of Theorem 5.6. The prin- 
cipal part of da/dL is 



4m / r<Ti 
'0 



ii ( ill I ii 



[)dx 



and 



m{r + 0'i) 



R{x) 

Therefore, it suffices to find an equivalent of 



Rdy. 



26*1 



R{x ~ L) 



R{y — V) dydx. 
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We replace (j)' by its expansion (6.2) which imphes that the expression that we wish 
to estimate is 



2e2^ Ind^- 



2e2^ ln(l 



1 + e^ 



■ dx 



C3 



■ dx 



1 



Each of the three integrals is equal, up to terms of order O(y^) to the following 
values: 

the computation is straightforward, and uses the change of variable 

1 + e^- = ^ 



1-2/ 



All these calculations lead to the following principal part: 



da _ 1 



3^2 27ji3 



64 



4 16 16 



+ 0{L-^^^). 



Finally, we find that 



Mn(V^) = -A^ + +{. + .')0{L-^'-). 



where x(/^) is given by 



(6.5) x{^^) 



r^2 



16 



1 



(2Ai2 - 15^3/8)^ 



4 16 

Thus, the skew stabilization takes place if 

> 

The critical value Vc for which the stabilization takes place is given by 

2x(m) 



(6.6) 



Vr ~ V- 



3L2 



We find the value given at proposition 5.8, when /i2 = 1 and the other /ij's vanish. 

When the coefficient m of —d^jdx^ is replaced by m + i^/eyLQ^ the skew stabi- 
lization also takes place under the same condition and for the same critical value 
of e, up to higher order terms. 

Therefore, we have proved the final result of this article: 



Proposition 6.1. In the general case i.e. for the equation 



ut 



■ iy/efio)uxx - {m + i^/e ^ii)u + (l + iVe, A^2)|upu - (l + «VeM3 



there exists a pulse-like solution for all s ma ll enough v, for all p > 0, and for all 
y < U' ; provided that x{fi) defined by (6.5) is strictly positive; in this case, the 
skew .stabilization takes place, at a critical parameter whose equivalent is given 
by (IJl). 
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7. Appendix. An implicit function Theorem with estimates 

The purpose of this appendix is to prove a version of the implicit function The- 
orem which is appropriate for singular perturbations. 

Lemma 7.1. Let X and Z he Banach spaces, and let f be a function from a 

neighborhood U of xq X to Z . Let zq = f{xo). Assume that A = Df{x{)) has a 
bounded inverse A~^ . Assume that the ball of radius p and of center xq is included 
in U. Let 

M= sup \\A-^D''f{xQ + 0\V 

l«l<p 

There exist constants a and K given by 

a = min(l, {2pM)-^), K=^ 

such that if \A~^zo\ < K, the equation 

fix) = 

possesses a unique solution in the ball {\x — xo\ < ap}; moreover, this solution 
satisfies 

\x - Xo\ < 2\A^^zo\ and \x - xq + A^^zq] < 2M\A~^zo\'^. 

Proof. Let 

i) = ^ - A-\f{xo + + (1 - t)A-'zo. 

It is equivalent to find a fixed point of F for t = 1 and to find a solution of f{x) = 0. 
Introducing F enables us to start from t = 0, where we have -^(0, 0) — 0. The 
implicit function Theorem in its classical formulation would give us the existence 
of an interval of t on which we can find a solution of ^ — t) = 0, but it would 
not ensure that we have a solution up to t = 1. The purpose of our estimates is to 
show that we can go as far as t = 1 . 

Let us determine a E (0, 1] such that ^ i—^ F{^, t) is a contraction of ratio 1/2 for 
ICI < ap and < t < 1: 

= ^'-^-A-' [ Df{xo + ^ + i(r - - dt 

Jo 



= i'-i-A-' 



Df{xo){e-0 



dt 



+ f D^f{xo + + 1{^' - 0)) + 1{^' - 0) (C' - ds 

Jo 

= - f f D^f{xo + s{^ + t{^' -i))){^ + t{^' -Odsdt. 
Jo Jo 

Therefore, we have the inequality 

m',t)-F{^,t)\<Map\^'-^\. 
Thus, the first condition that we wish to impose is 

Map < ^. 
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The second step is to impose a condition on a and on |i?~^2;o| such that ^ 
F{^, t) maps the ball of center xo and radius ap into itself, for t <1. We have 



Therefore, it is enough to require that 



Df{xo)^+ f D^f{xo + tOe^{l-t)dt 
Jo 



\A-^zo\ < 



Zap 



Thanks to the strict contraction Theorem, there exists for t € [0, 1] and for 
I a; I < ap a unique solution of 

Let us denote this solution by ^ = g{t). We may estimate g{t) as follows: we have 
the inequality 

\g{t)\ < \F{g{t),t) - F{0,t)\ + \F{0,t)\ < ^\g{t)\ +t\A-^zo\. 

Therefore, 

\g{t)\ < 2t\A-^zo\. 
But we can obtain a much better estimate, since 
g{t) + tA-^zo = F{g{t), t) + tA'^zo 

= g{t) - A-^f{xo + g{t)) + (1 - t)A-^zo + tA-^zo 
= g{t) - A-^za - A-^Df{xo)g{t) 

- [ D^f{xo + sg{t))g{tf^{l-s)ds + A-'zo 



I 

Jo 

7" 

Jo 



Therefore, 



D'f{xo + sg{t))g{tf^l-s)ds. 



M. 



\g{t)+tA-'zo\ < -j\9{t)f < 2M|A-i^o|', 



which concludes the proof. 



□ 
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